Astronomy & Astrophysics manuscript no. martins 


© ESO 2008 


February 5, 2008 





o 
o 

(N 



(N 
CO 

o 

O 

S3 

Oh 
6 

03 



X 



Stellar and wind properties of massive stars in the central parsec 

of the Galaxy 

F. Martins 1 , R. Genzel 1 - 2 , D.J. Hillier 3 , F. Eisenhauer 1 , T. Paumard 1 ' 4 , S. Gillessen 1 , T. Ott 1 , and S. Trippe 1 

1 Max-Planck Institiit fur extraterrestrische Physik, Postfach-1312, D-85741, Garching, Germany 
e-mail: martins@mpe.mpg . de 

2 Department of Physics, University of California, Berkeley, CA 94720, USA 

3 Department of Physics and Astronomy, University of Pittsburgh, 3941 O'Hara St., Pittsburgh, PA 15260, USA 

4 LESIA, Observatoire de Paris, 5 Place Jules Janssen, F-92195, Meudon, France 

Received 3 November 2006 / Accepted 1 March 2007 

ABSTRACT 

Context. How star formation proceeds in the Galactic Center is a debated question. Addressing this question will help us understand 
the origin of the cluster of massive stars near the supermassive black hole, and more generally starburst phenomena in galactic nuclei. 
In that context, it is crucial to know the properties of young massive stars in the central parsec of the Galaxy. 

Aims. The main goal of this study is to derive the stellar and wind properties of the massive stars orbiting the supermassive black hole 
SgrA* in two counter-rotating disks. 

Methods. We use non-LTE atmosphere models including winds and line-blanketing to reproduce H and K band spectra of these stars 
obtained with SINFONI on the ESO/VLT. 

Results. The GC massive stars appear to be relatively similar to other Galactic stars. The currently known population of massive stars 
emit a total 6.0 x 10 50 s _1 (resp. 2.3 x 10 49 s _1 ) H (resp. Hei) ionising photons. This is sufficient to produce the observed nebular 
emission and implies that, in contrast to previous claims, no peculiar stellar evolution is required in the Galactic Center. We find that 
most of the Ofpe/WN9 stars are less chemically evolved than initially thought. The properties of several WN8 stars are given, as 
well as two WN/C stars confirmed quantitatively to be stars in transition between the WN and WC phase. We propose the sequence 
(Ofpe/WN9 ^ LBV) — > WN8 — > WN/C for most of the observed GC stars. Quantitative comparison with stellar evolutionary tracks 
including rotation favour high mass loss rates in the Wolf-Rayet phase in these models. In the OB phase, these tracks nicely reproduce 
the average properties of bright supergiants in the Galactic Center. 

Key words. Stars: early type - Stars: Wolf-Rayet - Stars: atmospheres - Stars: fundamental parameters - Stars: winds, outflows - 
Galaxy: center 



5r i 1 ■ Introduction Galactic Center could be a peculiar environment for the forma- 
tion of massive stars. 

The center of our Galaxy is a unique environment to study mas- ,. . , , ■ 

sive stars. It harbors three of the most massive clusters of the As a P ecullar environment, the central cluster is especially 

Galaxy - the Arches, Quintuplet and central clusters. Heavily interesting since it harbors the supermassive black hole SgrA . 

extincted and only accessible at infrared (and longer) wave- J^^V^ofJ? obser , y, ed in mis region were the so-calle d 

lengths or in X-rays, each of these clusters has a population "f. ™* ™ Sl6 stars jForrest et alj B | Allen et alj llggj. 

of more than a hundred massive stars. Even more interesting is Their K " band s P ectmm showed stron 8 emlsslon llnes ~ ln P^ 10 " 

the difference in their ages: 2.5, 4 and 6 Myrs for the Arches, ul f * e H ^ ' 2.058 feature - and they were immediately clas- 

the Quintuplet and the cent ral cluster respectively dFiger et all srfied as Ofpe/ WN9 a class ofgy^lved massive stars Further 

1999, 2002; PaumaHUD [2006). Such a spread implies the observations by | Krabbe et al. | il99l|) revealed a few additional 

presence of different types of massive objects, naturally sam- He i emission line stars. They argued that this population most 



pling stellar evolution in the upper HR diagram. The youth of probably resulted from a burst of star formation a few Myrs 
these clusters, together with their total mass in excess of 10< M , a f° „S^,S^SE M 5 ; i 



partly explains the large number of massive stars in the Galactic atmosphere models by | Naiarro et al.| (EM EM) established 

Center. But another reason may be the top-heavy mass function: * at * ese oh ^ s were P ost :™ sequence massive stars Their 

Stolteltd] 42002) for the Arches and Paumard et al. ^M) for UV flux was able to . lomse th f ISM and P roduce * e nebular Br y 

the central cluster have shown that the slope T of the present-day emi ssion, but was far too soft to explain the nebular He i emis- 

mass function was shallower than the standard Salpeter value ^ on - Tbe P resence of a P<>P^on of hotter objects not accessi- 

(-0.8 instead of -2.35). This may be due to mass segregation ble to observations was thus inferred. 

or to a true feature of the initial mass function. In that case the This popu l ation was unraveled in the last years. 

Pau mard et alj d200ll) found additional stars with broader 

Send offprint requests to: F. Martins lines than the initial "He i" stars and typical of Wolf-Rayet stars. 
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iGenzel et al.l ((2003) deduced the presence of even more massive 
stars from the absence of CO absorption bands. The advent 
of adaptive-optics ass isted integral field sp ectroscopy lead to 
a new breakthrough: Paum ard et alj d2006l) spectroscopically 
identified nearly a hundred massive stars, including various 
types of Wolf-Rayet stars, O and B supergiants, and even dwarfs 
(see a lso Horrobin et al. 2004; Paumard et al. 2004). Previously, 
iGhez et al.l (120031) and Eisenhaue r et al.l d2005l) had shown that 
the group of stars located in the central arcsecond of the Galaxy 
and orbiting very close to the black hole (the so-called 'S-stars') 
was composed of early to late B dwarfs. 

This latter group of stars is at the heart of an issue usually 
referred to as the "paradox of youth" of the Galactic center: 
how could star formation (traced by the presence of young mas- 
sive stars) have happened so close to the black hole, where the 
tidal forces sh ould prevent any molecular cloud from collapsing 
(Morrisl lT993l) ? Two scenarios are invoked to solve this prob- 
lem. In the "in-situ formation" scenario, a disk forms around 
the supermassive black hole with a density large enough to 
be self-gravitating so that tidal f orces do not perturb the col- 
lapse of molecular materia l (e.g. iLevin &Beloborodovll2003t 
iNavakshin & Cuadrall2005l) . In the alternate "in-spiraling clus- 
ter" scenario, young massive stars are born in a dense cluster sev- 
eral tens of parsec away from SgrA* and its hostile environment, 
and are subsequently brought to the central regio n by spiral-in of 
the parent cluster du e to dynamical friction (e.g. Gerhard 2001; 
Kim & Morris 2003). The former picture (star formation in a 
self-gravitating disk) appears more attracti ve since it account s 
for a larger number of observational facts (Paum ard et al.l l2006). 
but the question is not completely settled. 
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Fig. 1. SINFONI K-band spectra of four post main sequence 
massive stars observed in the GC. The main lines are marked. 
Note that all lines are not present simultaneously in one spec- 
trum: depending on the stellar and wind parameters (in partic- 
ular T e ff and abundances), only a subset of the marked lines is 
observed in a given star. 



In that context, it is especially important to get a better 
knowledge of the physical properties of the massive stars in the 
central parsec of the Galaxy. In parallel to the development of 
powerful infrared observational techniques, reliable atmosphere 
models have become available in recent years. This is due to a 
huge effort from different groups to include thousands of metal- 
lic lines in such models ("line-blanketing"). The use of such 
models has quantitatively changed our knowledge of the stellar 
and wind properties of massi ve stars. To name a few, their effec- 
tive temperatures are lower dMartins et alj|2002t ICrowther et al.l 
120021: iRepolust et al.ll2004l) and their winds are highly c l umpe d 
dffillier et al.l 120031; feouret et all 120051; iFullerton et ail [2006). 
In addition to these significant developments, special atten- 
tion was given to the infrared range. Several studies have been 
undertaken to compare the parameters derived from IR di- 
agnostics with those obtained from classical opti cal features 
(Boh annan & Crowfherlll999l: IRepolust et ai]|2005l) . It appears 
that the differences are usually small, and certainly within the 
uncertainties on the derived parameters. Then, both new obser- 
vational data and better atmosphere models are now available. 
A quantitative investigation of the stellar and wind properties of 
massive stars in the central parsec of the Galaxy, in the context 
of the "paradox of youth", is thus possible. 

We present such a study in this paper. Sect. [2] describes the 
observational data. In Sect. [3] we present our atmosphere models 
and explain our method to derive stellar and wind properties. 
This is followed by a detailed analysis of individual stars (Sect. 
@). We then discuss the results with special emphasis on ionising 
fluxes (Sect. |5), stellar evolution (Sect. |6j, metallicity (Sect. |7), 
stellar winds (Sect. [8) and chemical composition of OB stars 
(Sect. [9). We finally summarize our findings in Sect.lTOl 



2. Observations 

The observations analysed here were conducted with the integral 
field spectrograp h SPIFFI/SINFONI on t he ESO/VLT Yepun 8 
meter telescope dEisenhauer et al.l 2003a b; [Bonnet et alj [2004) 
as part of the MPE-Garching GTO program "Galactic Center". 
A first mosaic of data cubes was obtained on Apr. 8 th 2003 us- 
ing the 250 mas scale in the K band mode, allowing a resolu- 
tion of ~ 4000. A second mosaic was observed in H+K band 
on the same scale on Apr. 9 th 2003. These observations were 
conducted in seeing limited mode since at that time the AO sys- 
tem MACAO was not yet coupled to SPIFF! A new SINFONI 
mosaic was obtained on Aug. 18-19 th 2004 in K band with the 
100 mas scale in adaptive optics mode. Finally, several fields of 
the region ~ 15 " north of SgrA* were observed on Mar. 16- 
17 th 200 5. Additional in f orma tion on all these data cubes can be 
found in iPaumard et al.l d2006l) . Two new fields were obtained 
in the configuration H+K/0. 1 mas scale on Apr. 20'' 1 2006 and 
Aug. 16"' 2006. 

To extract the spectra out of these cubes, we defined "source" 
pi xels showing spectral s ignatures of the massive stars identified 
by IPaumard et alj |2006) from which we removed neighboring 
"continuum" sources to correct for the local background. In the 
definition of the "continuum" pixels in crowded regions, we paid 
special attention not to include pixels contaminated by neigh- 
boring stellar sources. We used several combinations of source- 
continuum pixels to check the reliability of the extracted spec- 
trum. We found that as long as contaminating sources are not 
included, the spectra can safely be extracted. Fig. Q] shows typi- 
cal spectra of various types of stars analysed here together with 
line identification. 
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3. Modeling 

3.1. Atmosphere models 

In order to derive the stellar and wind properties of the massive 
stars in the central parsec of the Galaxy, we have used state of 
the art atmosphere mo dels computed with the code CMFGEN 
dHillier & MilleJ [T998). Such models include the main ingredi- 
ents necessary to produce realistic atmospheric structures and 
emergent spectra, namely a non-LTE treatment, winds and line- 
blanketing. The latter has been included only recently since in 
combination with the two former ingredients it leads to compli- 
cated and numerically demanding simulations. But the effects of 
line-blanketing are qualitatively and quantitatively very impor- 
tant. Most stellar parameters have to be revised when determined 
with line-blanketed models due to the strong modification of the 
radiative transfer caused by additional opacities from metals. 

CMFGEN computations proceed in two steps: first the atmo- 
spheric structure and radiative field are computed in an iterative 
process; second, a formal solution of the radiative transfer equa- 
tion with opacities given by the first step is determined, lead- 
ing to the detailed emerge nt spectrum. Details about the code 
CMFGEN are discussed bv lHillier & Miller! (I1998I) and here we 
only recall the main characteristics: 

- non-LTE treatment: populations of individual energy lev- 
els are computed through the resolution of statistical equi- 
librium equations including radiative and collisional pro- 
cesses. Non-LTE is especially important for infrared studies 
since in this range, stimulated radiative proce sses are greatly 
enhan ced compared to shorter wavelengths dLenorzer et al.l 
2004) and thus have a strong impact on level populations. 

- winds: a spherical geometry is adopted to fully take into 
account the atmospheric extension due to winds. Velocity 
gradients in the accelerating wind are also included in the 
radiative transfer problem. An important point to be noted 
here is the fact that a velocity law (equivalent to a density 
law through the equation of mass conservation) has to be 
adopted since radiative acceleration is not used to compute 
the hydrostatic structure of the atmosphere. In practice, two 
approaches are used: for OB stars, a photospheric structure 
comp uted with another atm osphere code (usually TLUSTY, 
see iHubenv & Land 1 1995b is smoothly connected to a so 
called "y8 velocity law (v = v«>(l - ^f-~f ) where v«, is the ter- 
minal velocity and /?* the stellar radius; for stars with denser 
winds, the atmosphere is usually optically thick so that the 
inner velocity structure is less crucial and a law of the form 

(voo-voxi-^y 3 

v = vo + j^— (1) 

1 + -^-e "elf 
Vcore 

is adopted (v core being the velocity at the bottom of the at- 
mosphere, vo the velocity at the expected photospheric ve- 
locity and /i e jf the density scale height of the photosphere)./? 
can in principle be derived from the s hape of spectral lines 
dPuls et alJ 119961; iMartins et all 120041) but here, it is only 
possible in a few cases. Hence, for m ost of the analysis , 
we simply adopt the standard value 1.0 (Na iarro et alJI 19941 
ICrowther et alj|2006l) . 

CMFGEN also allows the treatment of non homogeneous 
atmosphere through the adoption of a clumping law of the 
form 

/ = /«,+(!- U)e~^ (2) 



where is the value of / at the top of the atmosphere and 
Vi n it is the velocity at which clumping appears. We adopted 
a typical value of /«, = 0.1 for the remaining stars, unless 
explicitly indicated. 

- line-blanketing: CMFGEN includes a direct treatment of 
metals. The main approximatiorQ (which can be easily 
dropped provided the computational resources are available) 
is the grouping of levels of similar e nergie s in "super-levels" 
as initially proposed by Anderson (1991). Through the res- 
olution of the statistical equilibrium equations, level popu- 
lation of metals are computed as for H and He. In the mod- 
els presented here, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, 
Fe, Ni are included for Wolf-Rayet stars and related objects, 
while C, N, O, Si, S and Fe are used for OB stars. This is 
due to the lower density winds of the latter which require a 
better spatial sampling. In that case, minor metals have to be 
dropped to keep the size of the models reasonable. Note that 
no super-levels were used for H, He i and He n so that lines 
from these elements do not suffer from this approximation. 
In WC stars, we also did not use super-levels for C in and 
C rv. Finally, we tested the effect of super-levels on the N m 
2.247, 2.251 pm lines in WN stars and found it was n egligi- 
ble. In the following, we refer to lGrevesse & Sauvall (1998) 
for the solar abundances. 

- temperature structure: the temperature structure in the atmo- 
sphere is set under the condition of radiative equilibrium. 

- microturbulent velocity: a microturbulent velocity (v tur b) 
must be provided in the computation of both the atmospheric 
structure and the detailed emergent spectrum. In the former 
part of the simulation, v tul -b is independent of the position 
in the atmosphere and we usually choose typical values of 
20 km s _1 for OB stars and 50 km s _1 for Wolf-Rayet stars. 
For the computation of the detailed emergent spectrum, we 
adopt v tu ,b — 10 km s -1 . Note that we have run test models 
with larger values of v tur b for the computation of the emer- 
gent spectrum, and found no difference since the lines are 
mainly shaped by the wind terminal velocity, much larger 
than v,^. 

3.2. Method 

The main diagnostics we used for the analysis were: 1) K band 
photometry and 2) normalized spectra. The number of diagnos- 
tics is relatively limited, which sometimes results in degenera- 
cies in the derived parameters, especially when T e ff cannot be 
accurately constrained. The procedure to derive the stellar and 
wind parameters is as follows: 

- Effective temperature: we relied on the ratio of lines from 
successive ionisation states of a given element: He for OB 
and WN stars, He and C for WN/C and WC stars. In practice, 
the following lines were used: He 1 2.1 12 pm, He 1 2.184 pm, 
Hen 2.037 pm, Hen 2.189 pm, Hen 2.346 pm, Civ 2.070 
pm, C iv 2.079 pm, C iv 2.084 pm, C in 2.325 pm. The He i 
line at 2.058 pm has been shown to gre atly depend on basi- 
cally any parameter (Naiarro et al Jl 1994b and, although it is 
frequently the strongest observed line, it is not used to de- 
rive T e ff. Unfortunately, two successive ionisation ratios of 
the same element are not always present, especially for the 
latest OB and Wolf-Rayet stars. In that case, one usually sees 
only He i lines so that only an upper limit on T e ff can be de- 
rived. 

1 Other approximations concern the treatment of line profiles, redis- 
tribution functions and microturbulence. 
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We also define here a temperature T+ such that 

L = AnRlTi = 4nR 2 2/3 T* B (3) 

R+ is the radius at which the Rosseland optical depth 
TRosseiand is equal to 20 (R2/3 corresponding to the radius 
where the TR OSSe i anc j = 2/3). This definition is useful in the 
case of evolved massive stars since it is not directly affected 
by the stellar wind, and since it allows a better comparison 
of stellar parameters to evolutionary models (see Sect. |7.2| >. 

- Luminosity: the main constraint on the luminosity comes 
from the absolute K magnitude (Mk) which is derived from 
the observed nix , the distance to the G alactic Center (7.62 
kpc according to lEisenhauer et al.ll2005l). and the extin ction 
taken from the recent work of ISchoedel & et all (120071) (see 
Table [TJ. Knowing the distance of the sources is a great 
advantage over many Galactic studies of massive stars. 
The main source of uncertainty in the luminosity is the 
extinction which is known to vary in the central parsec. 
A given value of Mk is obtained for a given flux in the K 
band which depends on 1) the luminosity and 2) on the 
wind density. Indeed, free-free emission in the atmosphere 
can produc e an excess of emission mimicking a larger 
luminosity dLamers & Cassineilil|1999l) . Hence, luminosity 
is derived in combination with wind parameters (M and Voo). 

- Mass loss rate: M is derived from the strength of emission 
lines formed by recombination. Their intensity depends on 
the wind density but also on the ionisation of the atmosphere 
(which is also partly controlled by the density) and on 
the abundances. The mass loss rate is thus derived in 
combination with effective temperature and abundances. 

- Terminal velocity: In stars showing blueshifted absorption 
profiles and P-Cygni lines, v M was estimated from the 
velocity shift of the bluest part of the absorption trough with 
respect to the rest-frame wavelength. For Wolf-Rayet stars 
with broad emission lines, half of the line width was chosen 
as typical of the terminal velocity. Values were then refined 
according to the quality of the fit. 

- Abundances: The relative H to He abundances can be 
derived from the intensity of various lines from these 
elements. As mentioned previously, such a determination 
cannot be separated from the estimate of M and T e ff so that 
a simultaneous determination of these parameters is done. In 
stars showing C and N lines (C iv 2.070 /im,Civ 2.079 fim, 
Civ 2.084 fim, Cm 2.325 fim, Nm 2.247 fim, Nm 2.251 
fim), constraints on carbon and nitrogen abundances can 
also be given. Note that we do not use the feature at 2.115 
fim for the N abundance deter mination: it is unclear whether 
N in is the only contributor dGeballe et alj 120061) and, in 
addition, the transition oscillator strength is uncertain. 

- Mass: The determination of masses for hot stars is a 
challenge. Here , we h ave used the M - L relation of 
Heger & Langei] ( 1 19961) to estimate the present mass of H 
free Wolf-Rayet stars. For other stars, we did not try to 
derive masses, the uncertainties being too large. 

The uncertainties in the derived parameters are the follow- 
ing: ± 3000 K for T eff (+ 6000 K for the Ofpe/WN9 stars), + 
0.2 in log 0.2 dex for M, 30 % for the abundances. They are 



not statistical errors (the estimate of such errors would imply the 
computation of a huge number of models to sample the param- 
eter space around the best fit solution). Instead, they reflect the 
range of values leading to acceptable fits of the observed spectra. 
Our errors also do not include any systematic contribution due to 
uncertainties in atomic data such as collisional and dielectronic 
recombination cross-sections. This has to be kept in mind when 
considering the values we quote. 

4. Analysis of individual stars 

In this section, we give the results of the detailed analysis of in- 
dividual stars. We mainly focus on Wolf-Rayet and the so -called 
"He i" stars since their spectra have large enough S/N ratios to 
allow quantitative spectroscopy. OB supergiants are studied by 
mean s of the average spectrum of 10 of them (Paumard et al. 
2006). The global methodology presented in Sect. I3.2l is not re- 
peated for each star: we only give specific comments when nec- 
essary. 

However, as a preamble, we would like to say a few words 
about the behavior of He i 2.058 fan. It is well known that this 
line is extremely sensitive to any detail of the modeling, and in 
particular to the amount of UV radiation. Indeed, the upper level 
of He i 2.058 fim is directly coupled to the ground state by a 
transition at 584 A. Dramatic improvement has been achieved in 
recent years in the prediction of this UV radiation, mainly due to 
the inclusion of line-blanketing in the models. However, as we 
will see in the following analysis, H ei 2.058 fim is still p oorly 
reproduced in several stars. Recently, Naia rro et al.1 J2006) have 
shown that the radiation at 584 A was partly controlled by the 
strength of two Fe iv lines: artificially changing the strength of 
these lines improved the fit of optical singlet Hei lines. We have 
tried the same kind of tests in the present study, but it turned out 
that the resulting He 1 2.058 fim line profiles were little changed. 
The reason is partly that in the Wolf-Rayet stars studied here, 
He i 2.058 fim is an emission line and is controlled b y recombi- 
nation processes, while in the O stars analysed by Naiarro et al. 
(HH), He i 2.058 fan is in absorption and depends much more 
on pure radiative transfer effects. Besides, the optical depth of 
the Hei 584 A line is large in Wolf-Rayet stars, partly control- 
ling the population of the He i 2.058 fim upper level. This does 
not mean that Hei 2.058 fim is not sensitive to the UV contin- 
uum: as we will discuss in Sect. [7] subtle blanketing effects can 
significantly change the appearance of He i 2.058 fim. Having 
said that, we now turn to the detailed study of individual GC 
stars. 

4.1. Ofpe/WN9 stars 

In this Section we present the analysis of five Ofpe/WN9 stars. 
We deliberately exclude the stars I RS 1 6S W and IRS 1 6 NE since 
they are binaries (or candidates, see lMartins et alj|2006l) . 

4.1.1. IRS34W 

IRS34W is an Ofpe/WN9 star and the fain test of the LBV can - 
didates iden t ified i n the Galactic Center (Pau mard et al.l l2004). 
Tripp e et"aT] (|2006) also showed that it had recently experienced 
photometric variability attributed to the formation of dust in ma- 
terial possibly ejected in a LBV-like outburst. 

Due to the absence of a strong T e ff indicator, we could find 
several solutions to the fit of the K band spectrum by varying the 
He content, the luminosity, and the mass loss rate for T e ff in the 
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Table 1. Observational properties of the stars an alyzed in this paper. Spectral types are from Paumard et al.l d2006l) . Values of 
extinction are taken from lSchoedel & et al.l d2007l) . A distance of 7.62 kpc is adopted (Eise nhauer et alj|2005l) . 
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-6.0+0.3 
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-4.4+0.3 
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12.8+0.1 
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-4.6+0.3 




16SE2 


WN5/6 


12.0+0.1 


2.6+0.2 


-5.0+0.3 


GCIRS 16SE2 


7W 


WC9 


13.1+0.1 


2.5+0.2 


-3.8+0.3 


GCIRS 7W 


7SE 


WC9 


13.0+0.1 


2.6+0.2 


-4.0+0.3 


GCIRS 7SE 


13E4 


WC9 


11.7+0.1 


2.8+0.3 


-5.5+0.4 


GCIRS 13E4 



range ~ 20000 — 33000 K. Fig. [2] shows one of the possible 
best fit models. For the acceptable effective temperatures, the 
He/H ratio is larger than solar, in the range 0.25 - 0.6, while the 
mass loss rate goes from 6 x 10~ 6 - 2 X 10" 5 M yr~' and the 
luminosity is found in between 3 and 6x 1O 5 L . It is important to 
note that these parameters are not independent and only certain 
combinations within these ranges are acceptable. For example, 
low T e ff imply low luminosities, large M and large He/H ratio. 

The morphology of Br y is well sampled in our observed 
spectrum. In particular, the blue side of the Br y emission line 
shows a "shoulder" and, at even shorter wavelengths, a small 
absorption dip. This results from the combination of two effects: 
first, Br r is a P-Cygni profile for which the blue absorption dip 
is partially filled by emission; second, several He lines with both 
emission and absorption profiles add to this Br r shape. The con- 
tribution of each element is shown in the insert in Fig. [2] The 
exact shape of this complex profile depends on the wind density 
in the line formation region, and consequently depends on both 
f3 and the filling factor / (in addition to the mass loss rate and 
the He/H ratio). In order to fit this profile, we had to choose ft in 
the range 2.0 - 4.0, i.e. larger than 1, the value commonly found 
in O supergiants. This is due to the fact that larger ft leads to nar- 
rowe r profiles with stron ger absorption dips and emission peaks 
(see iMartins et al.l 12004). Similarly, the best fits were obtained 
with unclumped models. This is at odds with the current knowl- 
edge that winds of massive stars are strongly inhomogeneous. 
However, since this complex line profile depends on several pa- 
rameters, we refrain from concluding that the wind of IRS34W 
(and the other similar Ofpe/WN9 stars) are homogeneous. More 
detailed investigations with high resolution spectroscopy should 
help to better resolve this line and improve our determination. 

The only feature not well reproduced by our models is the 
emission at 2.112-2.115 /mi: it is wider in the observed spec- 
trum than in our models. The main reason is that this emission 
is a blend of several lines, the identification of which is still 
under debate. The concensus is that both Hei 2.112 fim (blue 
part, through P-Cygni profiles) and N m 2. 1 15 fmi (red part) con- 
tribute to the emission. Clearly, the He i emission is present in 
our models but we are not able to reproduce the N m emission. 
Actually, we are not convinced that this emission is really due 
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Fig. 2. Best fit (red dot-dashed line) of the observed K band spec- 
trum of IRS34W (Ofpe/WN9, black solid line). The insert is a 
zoom on the Br y line, showing the contribution of H (blue dashed 
line) and He i (green dotted line) to the total synthetic profile (red 
dot-dashed line). 



to N in sin ce we do not detect the lines N m 2.247 fim and N m 
2.251 fim. iGeballe et al.l d2006l) recently claimed that O m could 
be responsible for that emission. Since the oscillator strength for 
this transition is very uncertain, we decided not to fit the red part 
of the emission complex at 2.1 12-2.1 15 //m. 

4.1.2. IRS16NW 

The spectrum of IRS16NW is similar to IRS34W with the ex- 
ception that the emission part of the 2.112 fim complex is weaker 



Table 2. Derived stellar and wind parameters. The typical errors are: ±3000 K on temperatures (+6000 K for Ofpe/WN9 stars), ±0.2 dex on log and logM, 100 km 
terminal velocities and ±30% on abundances (except special cases; see comments on individual stars). 



Star 


ST 


T, 


T c ff 




R, 


R2/3 


M K 


log M 


foo 




H/He 


C/He 


X(N) 


M 


log Qh 


log 0HeI 






[K] 


[K] 




[Ro] 


[Ro] 




[M yr->] 




[kms" 1 ] 


# 


# 




[Mo] 


[s" 1 ] 


[s- 1 ] 


34W 


Ofpe/WN9 


23000 


19500 


5.5 


35.9 


49.3 


-6.11 


-4.88 


1.0 


650 


4.0 








48.32 


47.64 


16NW 


Ofpe/WN9 


20000 


17500 


5.9 


59.1 


75.8 


-6.80 


-4.95 


1.0 


600 


5.0 








48.04 


47.35 


16C 


Ofpe/WN9 


21500 


19500 


5.9 


63.9 


79.5 


-7.05 


-4.65 


1.0 


650 


2.5 








48.76 


47.60 


33E 


Ofpe/WN9 


20000 


18000 


5.75 


63.9 


75.9 


-6.85 


-4.80 


1.0 


450 


4.0 








48.26 


47.60 


AF 


Ofpe/WN9 


23000 


21000 


5.3 


28.1 


34.5 


-5.77 


-4.75 


0.1 


700 


2.0 








48.06 


47.82 


15NE 


WN8 


34500 


33000 


5.25 


11.9 


13.0 


-4.65 


-4.70 


0.1 


800 


0.0 


<1 10- 4 


0.0137 


12.6 


48.56 


47.47 


AFNW 


WN8 


37000 


33000 


5.5 


13.9 


17.5 


-5.19 


-4.50 


0.1 


800 


0.1 


<1 10- 4 


0.0326 


18.6 


49.14 


47.60 


9W 


WN8 


40500 


32000 


5.4 


10.2 


16.4 


-5.20 


-4.35 


0.1 


1100 


0.1 


<5 10~ 5 


0.0133 


15.8 


49.11 


47.79 


7E2 


WN8 


37500 


34500 


5.2 


9.5 


11.1 


-4.25 


-4.80 


0.1 


900 


0.0 


<8 10~ 5 


0.0137 


11.7 


48.88 


47.70 


13E2 


WN8 


29000 


29000 


6.1 


44.7 


45.0 


-6.55 


-4.35 


0.1 


750 


0.1 


<3 10- 4 


0.0167 


82.5 


49.49 


47.61 


7SW 


WN8/WC9 


34500 


33000 


5.55 


16.8 


18.0 


-5.19 


-4.70 


0.1 


900 


0.0 


0.005 


0.0135 


18.6 


49.15 


47.66 


15SW 


WN8/WC9 


39000 


35000 


5.1 


7.9 


9.7 


-4.36 


-4.80 


0.1 


900 


0.0 


0.013 


0.0229 


10.3 


48.77 


47.83 


AFNWNW 


WN7 


36500 


28500 


5.25 


10.7 


17.2 


-4.59 


-3.95 


1.0 


1800 


0.1 


<1 10- 4 




12.6 


48.87 


47.77 


34NW 


WN7 


34000 


33000 


5.6 


18.2 


19.3 


-4.71 


-5.30 


0.1 


750 


1.0 


1.5 10~ 4 


0.0069 




49.25 


48.13 


16SE2 


WN5/6 


53000 


41000 


5.45 


6.4 


10.5 


-5.03 


-4.15 


0.1 


2500 


0.0 


<1 10~ 4 




17.2 


49.24 


48.51 


7W 


WC9 


47500 


39000 


5.1 


5.3 


7.8 


-3.79 


-5.0 


0.1 


1000 


0.0 


0.06 




10.3 


48.86 


47.80 


7SE 


WC9 


44500 


36500 


5.15 


6.4 


9.5 


-4.03 


-4.90 


0.1 


1000 


0.0 


0.04 




11.0 


48.88 


47.82 


13E4 


WC9 


42500 


37500 


5.8 


14.7 


18.7 


-5.44 


-4.30 


1.0 


2200 


0.0 


0.02 




45.0 


49.56 


48.47 
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Fig. 3. Best fit (red dot-dashed line) of the observed K band spec- Fig. 4. Best fit (red dot-dashed line) of the observed K band spec- 
trum of IRS 16NW (Ofpe/WN9, black solid line). trum of IRS 16C (Ofpe/WN9, black solid line). 



and the He lines on the blue wing of Br r are mainly in absorp- 
tion. However, IRS 16NW is 2.5 mag brighter, leading to a higher 
luminosity. The best fit model of its K band spectrum is shown 
in Fig. [3] For this star too, the only problem is the emission at 
2.112 /mi which is absent in our model. As for IRS34W, large /3 
and unclumped models gave the best fits. 

Compared to the analysis of Naia rro et alj {1997), we find a 
larger T e ff, although the range of values for which a fit can be 
achieved encompasses their value. The luminosity is lower by 
a factor 2.4, while M is lower by a factor of 4.7. We also find 
a ter minal velocity smalle r by 150 km s . The main difference 
with iNaiarro et al.l d 19971) is the He content that we find larger 
than solar but smaller than the H content. In our case, this in- 
dicates that IRS16NW has probably only recently left the main 
sequence and is in an early evolved status. 



4.1.3. IRS16C 

IRS16C is the brightest star of our sample. It is one of the 
Ofpe/WN9 stars first discovered in the GC. As for the other stars 
of the same spectral type, its T e ff is poorly constrained. Fig. [4] 
shows one of the best fit models. A large value of /3 as well as 
an unclumped wind are favoured as for IRS34W. The feature at 
2.115 fim is not reproduced in our model, as well as the one at 
2.100 /mi. These lines could be from Nm, but the absence of 
N in 2.247 /mi and N m 2.25 1 /mi emission weakens this possi- 
bility. We detect a Mgii emission at 2.138 //m. This emission is 
reproduced by our model for a twice solar Mg content. However, 
this value should be interpreted with care given the uncertainty 
in T e ff and the weakness of the line. 

As for IRS16NW, we find a lower He content and a low er 
mass loss rate compared to the results of Nai arro et al.l d!997l) . 




Fig. 5. Best fit (red dot-dashed line) of the observed K band spec- 
trum of IRS33E (Ofpe/WN9, black solid line). 



4.1.4. IRS33E 

The spectrum of IRS33E is similar to IRS34W, IRS16NW and 
IRS16C, and so are the derived properties. The best fit model is 
shown in Fig. [5] As for IRS16C, the weak emission lines at 2.100 
/mi and 2.1 15 //m are not reproduced. Note that the structure of 
the Br y emission is real and is observed at different epochs. 
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Fig. 6. Possible fit (red dot-dashed line) of the observed H+K 
band spectrum of the AF star (Ofpe/WN9, black solid line). 
Other combinations of T e ff , He/H and M lead to similar fits. The 
observed spectrum around 1.85 //m is noisy and should be disre- 
garded. See text for discussion. 



4.1.5. AF star 

The AF star is one of the first GC stars discove red and anal- 
ysed. B ased on its SINFO NI K band spectrum. | Fprrest et alj 
d 19871) . lAllen etalJ dl990b and iNaiarro et al] dl994j) classified 
it as a Ofpe/WN9 star. However, compared to the other stars 
of the same type presented above, AF has much broader lines. 
This indicates a stronger wind and possibly a more advanced 
evolutionary state (see also Sect. |6.2| i. The best fit to the H+K 
band spectrum is shown in Fig. [6] The effective temperature 
is poorly constrained due to the absence of Hen lines, as for 
the other Ofpe/WN9 stars. This implies a degeneracy between 
T eff and He/H: models with 18000/T < T eS < 30000/T and 
0.5 < He/H < 5.0 gave reasonable fits (lower He/H being 
required at larger T e g). Such behavior was already noted by 
iNajarro et al] (ll994). We also found that M ~ 10 4 M Q yr~' were 
necessary to fit emission at low T e ff, while values of the order 
1.5 x 10~ 4 M yr~' were sufficient at high T e ff. The luminosity is 
in the range 1 - 2 x 1O 5 L . The values derived bv lNaiarro et alj 
(1994) are consistent with our cool/He-rich/large M models. The 
analysis of this star illustrates perfectly the degeneracy one has 
to face when T e ff is poorly constrained. In contrast to the other 
Ofpe/WN9 stars analysed before, the shape of the Br r complex 
is dominated by the wind, so that the contribution of the He lines 
blueward of Br r is not resolved. Hence , the He content is p oorly 
constrained. In the comparison to the N aiarro et alj d 19941) re- 
sults, it is also important to note that a different extinction was 
used. We adopted Ak = 2.2 while Najarro et al. had Ak — 3.0. 
This influences the results, in particular the luminosity and the 
mass loss rates (through their influence on Mk, see Sect. 13.21 ). 



4.2. WN stars 

In the following, we present the results of the detailed study of 5 
WN8, 3 WN7, 1 WN5/6 and 2 WN8/WC9 stars.For all the WN8 
stars, the H content is relatively low, but we cannot discriminate 
between H free stars and stars with X(H) of a few percent. We 
thus quote a H mass fraction < 0. 1 . 

4.2.1. IRS15NE 

Fig. shows our best fit spectrum of the WN8 star IRS15NE. 
The Hen lines at 2.037, 2.189 and 2.346 /jm are well reproduced 
and allow a good estimate of the effective temperature. The fit of 
the N m 2.247, 2.251 ymi features is good, providing an accurate 
N abundance determination. The main discrepancy concerns the 
He i 2.058 fim line which is too weak. Reducing T e (j improves 
the fit but weakens the He n lines. The Si iv feature at 2.427 /mi 
is also slightly too weak if a solar abundance is used for Si. Fig. 
|7]shows the effect of increasing the Si content by a factor of 2.5 
and 7. Interestingly, the fit of the SilV line improves, as well as 
the fit of He i 2.058 fiml However, several weak Si lines appear 
around 1.98 and 2.08 ym. We do not detect these lines. Hence, 
we cannot safely conclude that a super-solar Si abundance is re- 
quired for IRS 15NE. The problem of Si iv 2.427 yum in the initial 
model may be due to incorrect atomic data for this line. In terms 
of atmospheric structure, changing the Si content translates into 
a slight variation of the temperature structure (T is reduced in 
the outer atmosphere) which is then responsible for the strength- 
ening of He 1 2.058 ymi (see also Sect. |7.21 >. This highlights once 
more the extreme sensitivity of He i 2.058 ym to the very details 

of the modeling. 

IRS15NE was previously studied by INaiarro etail (1 19971) 
who found a lower T e ff, a larger luminosity and a larger mass 
loss rate. The terminal velocity and He content were similar to 
the present value. 

4.2.2. AFNW 

AFNW is located North West of the AF star and was assigned 
a spectral type WN8 by Pau mard et al.l d2006l) . The presence of 
both He i and He n lines allows a rather robust T e ff determination 
(see Fig. |S). Only He i 2.058 ym is too weak and He i 1 .700 ym 
too strong in our model. 

4.2.3. IRS9W 

IRS9W is the WN8 star of our sample with the cleanest spectrum 
and the strongest Hen lines. Fig.|9]shows that our best model is 
able to perfectly reproduce most of the features, with the no- 
table exception of He i 2.058 ym. Given the quality of the ob- 
served spectrum and the presence of several Hei and He n lines 
very well reproduced by our model, we think the effective tem- 
perature is well constrained (see in particular the ratio of He n 
2.189 ym to He i 2.184 ym). Decreasing T e ff to strengthen Hei 
2.058 ym leads to weaker He n 2.037, 2.189 and 2.346 ym lines. 
Modifying the luminosity and mass loss rate so that they still 
lead to a good match of the observed K band magnitude and of 
the strength of emission lines does not lead to any improvement 
as far as Hei 2.058 ym is concerned. Interestingly, the situation 
got better when we tried to increase the global metallicity to a 
value of 2 times Z . In that case, we could get a strong He 1 2.058 
ym line without degrading the fit of the other diagnostics. This 
points once again to the extreme sensitivity of Hei 2.058 ym 
to UV opacities, a larger metallicity corresponding to a softer 
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Fig. 7. Best fit (red dot-dashed line) of the observed K band spec- 
trum of IRS15NE (WN8, black solid line). The different broken 
lines indicates models with similar parameters except the Si con- 
tent. When it increases, the fits of Si iv 2.427 fim and He i 2.058 
fan are improved. See text for discussion. 



Fig. 9. Best fit of the observed K band spectrum of IRS9W 
(WN8, black solid line). The blue dotted line is a model with 
Z = Z while the red dot-dashed line is for Z — 2 x Z Q . Note 
the change of He 1 2.058 fim with metallicity, tracing its extreme 
sensitivity to UV radiation. 




Fig. 8. Best fit (red dot-dashed line) of the observed H band (left) 
and K band (right) spectrum of AFNW (WN8). 



4.2.4. IRS7E2 

The best fit model to the K band spectrum of IRS7E2 is shown 
in Fig. [10] All lines are well reproduced, except Hei 2.058 fim. 
The excellent fit of the Nm 2.247, 2.251 fim line allows a good 
N abundance determination. T e fr is also well constrained since 
Hen 2.189 fim, Hen 2.037 fim and Hen 2.346 fim are clearly 
detected. As for IRS9W, we have tried to increase the metallicity 
to improve the fit of He i 2.058 fim, but this time, the line barely 
reacted and remained too weak. 



4.2.5. IRS13E2 

IRS13E2 is the brightest member of the IRS13E cluster and is 
classified as WN8. Our best fit model is shown in Fig. [TTJ An 
effective temperature of 29000 K was required to fit the He spec- 
trum. It is the most luminous WN8 star of our sample However, 
the presence of dust in the IRS13E cluster may hamper our de- 
termination. We will argue in Sect. l4.3l and l6.5l that we probably 
der ive only an upper lim it on the luminosity. 

iNajarro et al.l (Il997l) already mentioned the stars of IRS13E 
in their study, but at that time they could not resolve the compo- 
nents and analysed the global spectrum. Hence, a direct compar- 
ison with their results is meaningless. 

4.2.6. IRS7SW 



UV radiation. A more detailed discussion of the behavior of He i 
2.058 fim as regards metallicity changes is given in Sect. I7.2I 



iPaumard et all (I2006I) classified IRS7SW as WN8, but in view 
of the present results (see below), we refine its identification 
to WN8/WC9. Our best fit model is presented in Fig. [T5] The 
preferred T e ff allows a reasonable fit of the Carbon lines and of 
He ii 2.189 fim, but seems a little too low to account for Hen 
2.037 fim and Hen 2.346 fim. However, increasing T e ff leads to 
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Fig. 11. Best fit (red dot-dashed line) to the observed K band 
spectrum of IRS13E2 (WN8, black solid line). 



a too strong Hen 2.189 yum line. The presence of Civ and Cm 
as well as Nm lines indicates that IRS7SW is likely a WN/WC 
star. This is confirmed by the abundance determination (see also 
Sect. l6.3l for a quantitative discussion). We note too unidentified 
lines: one feature at 2.140 yum, and another one at 2.224 jjm. 
The former cannot be attributed to Mgn 2.138 fim since T e fr is 
too high. Besides, Mgn 2.138 ymi is a doublet while we clearly 
see only one component. The line at 2.224 yum cannot be due 
to Nai since again T e jf is too large. Interestingly, these lines are 



also found in most of the WC9 stars of our sample, as well as in 
the WN8/WC9 star IRS15SW (see next Section). We conclude 
that they are typical of C-rich stars. Note that these li nes are also 
presen t and not identified in the WC9 stars of the iFiger et ail 
( 1997) sample (see their Fig. 9). 

4.2.7. IRS15SW 

IRS15SW is a late WN star showing C lines in it s K band spec- 
trum s o that it was classified as WN8/WC9 by Paumar d et alj 
(2006). Most of the lines are reproduced by our best model (Fig. 
[T3l . A notable exception is the He i/He n complex around 2. 185 
/im. The emission is stronger than our model. Changing T e g does 
not help since we fit either the blue or red side of the emission 
but never the whole complex. Besides, T e ff is relatively well con- 
strained by the other He n lines. We also do not perfectly fit the 
blue absorption dip of Hei 2.058 /im. This may require a larger 
Voo, but in that case the other emission lines are too broad. The C 
lines are perfectly matched, allowing a reliable abundance (and 
T e ff) determination. IRS15SW will be further discussed in Sect. 
IO] 

IRS15SW was studied by iNajarro etall dl997l) . We find a 
much larger T e fr, essentially because we can rely on several He n 
lines and on the C iv/C m ratio (note the good fit of C lines in Fig. 
[T31 , We also find a much lower luminosity and a lower clumping 
corrected mass loss rate. The terminal velocity and H content are 
however similar (see Fig. [13] for a complete view of the He and 
H spectrum). 

4.2.8. AFNWNW 

AFNWNW is a WN7 star. Most of its lines are reasonably well 
reproduced by our best fit model (see Fig. \TM . The main prob- 
lems are the too weak Hei 2.058 /urn line and the too narrow 
Hei/Hen complex at 2.18-2.19 //m. Note however that the S/N 
ratio is rather low, preventing an accurate determination of the 
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Fig. 13. Best fit (red dot-dashed line) to the observed K band Fig. 15. Best fit (red dot-dashed line) to the observed K band 
spectrum of IRS15SW (WN8/WC9, black solid line). spectrum of IRS34NW (WN7, black solid line). 




Fig. 14. Best fit (red dot-dashed line) to the observed H+K band 
spectrum of AFNWNW (WN7, black solid line). The region 
around P a was cut since it is not reliable. The observed spec- 
trum was also smoothed for clarity. 



physical parameters. The absence of C iv emission indicates an 
upper limit of 10" 4 for C/He. Note that our best fit model is 
unclumped. We cannot derive the clumping factor from the ob- 
served spectrum due to the poor S/N ratio. 



4.2.9. IRS34NW 

IRS34NW is a WN7 star. It has narrower and weaker lines than 
AFNWNW, reflecting its weaker wind. The C iv lines indicate a 
slightly sub-solar C abundance (see best fit in Fig. [15]). Together 
with the N content derived from N m lines, this reveals an early 
stage of CNO processing (compared to AFNWNW which has no 
detectable Civ lines). IRS34NW still has a significant amount 
of Hydrogen. We conclude that IRS34NW is less evolved than 
AFNWNW in spite of a similar spectral type. This likely reflects 
different initial masses. 



4.2.10. IRS16SE2 

IRS16SE2 is the ear liest WN star of our sample (WN5/6, see 
Paumard et al. 2006). Consequently, it also has the highest ef- 
fective temperature (41000 K) which is quite well constrained by 
the presence of both He i and He n lines in the K band spectrum. 
Helium is indeed doubly ionised in most of the atmosphere, ex- 
cept in the very outer part where it recombines, leading to the 
He i absorption trough near 2.044 fim. This feature being due 
to a blueshifted Hei 2.058 pmi absorption, we have a good esti- 
mate of the wind terminal velocity (2500 km s _1 ). The absence 
of C iv emission around 2.08 fim - expected for such a large T e jf 
- sets an upper limit on the Carbon abundance (C/He < 10~ 4 by 
number), indicating CNO processing. 

ICrowther & Smithl d 19961) studied two WN6 stars with K 
band spectra very similar to IRS16SE2. Their results are in ex- 
cellent agreement with the present ones: ~ 55000/f, log ^ ~ 

5.4 and log 4= ~ -3.9. Since ICrowther & Smithl dl996h found 

no difference between their IR analysis and optical results, we 
are confident that the parameters we derive for IRS16SE2 for 
pure IR diagnostics are reliable. 
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Fig. 16. Best fit (red dot-dashed line) of the observed K band 
spectrum of IRS 16SE2 (WN5/6, black line). 

4.3. WC9 stars 

Here, the stellar and wind parameters of three WC9 stars are 
derived. One important word of caution is necessary though . 
WC9 stars are often associated with dust (Willia ms et alj fl987). 
The origin of this dust is not completely understood, but wind- 
wind interaction in binary systems is the favoured mechanism. 
Observation of dust spirals (also named 'pin whe els') around 
sever al WC9 stars strongly support this scenario (Tuthi lTet al.l 
Il999l) . Recent observations of the 'Cocoon stars' after which 
the Quintuplet clu ster was named showed such pinwheels 
(Tuthi ll et al.ll2006l) . 

The presence of dust in WC9 stars complicates the analysis 
of their IR spectra since it produces an additional emission which 
adds to the stellar-)- wind continuum. In practice, if WC9 stars are 
analysed under the assumption that they are dust-free, their con- 
tinuum is over-estimated. Consequently, when normalizing their 
spectra, lines appear weaker. This implies under-estimates of the 
mass loss rates and abundances. In addition, the derived lumi- 
nosity is over-estimated since the total continuum is composed 
of both the stellar+wind continuum and the dust emission. Note 
however that T e ff estimates are less affected, since the ratio of 
lines is only weakly affected by the presence of dust. 

In the following, we discuss for each star the observational 
evidence for dust and the reliability of the derived parameters. 

4.3.1. IRS7W 

L-band observations of IRS7W were recently performed by 
iMoultaka et al. (2005)0. Inspection of their Table 1 shows that 
the colors of IRS7W are consistent with the extinction law of 
the GC, at least in the HKL bands. A possible excess emission is 
only seen in the M band due to a red L-M color. Consequently, 
we think that our modelling, restricted to the H+K band spec- 

2 IRS7W is their WR2 star, and not WR1 as they claim. Their WR1 
star is IRS7SW. 
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Fig. 17. Best fit (red dot-dashed line) of the observed H and K 
band spectra of IRS7W (WC9, black solid line). 



trum, is not hampered by any dust emission. The derived param- 
eters can be trusted within their error bars. 

Our best fit is shown in Fig. [T7] The effective temperature 
is relatively well constrained by the presence of several Hen 
lines as well as C iv and C m features. The main discrepancy is 
once again the too weak He 1 2.058 [im in our model. In addition, 
we note the two unidentified lines at 2.140 and 2.224 fim as in 
WN/C stars. 

4.3.2. IRS7SE 

IRS7SE is a WC9 star very similar to IRS7W. Unfortunately, we 
do not have any information on its photometry (except for the K- 
band). We are thus not able to check the possible contamination 
by dust. Adopting a conservative approach, we consider that our 
results are only limits (lower for M and the C abundance, up- 
per for log j^). Fig. rT~8l shows a fit of similar quality compared 
to IRS7W, with the same caveats (Hei 2.058 fan, unidentified 
lines). 

4.3.3. IRS13E4 

Mail lard et al.l (l2004l) presented a detailed investigation of the 
stellar content of IRS13E. Using HKL photometry, they report 
the discovery of several very red sources in addition to the 
bright components IRS13E2 and IRS13E4. They interpret these 
sources as dusty Wolf-Rayet stars. Unfortunately, IRS3E4 is not 
detected in their L band images, so that the presence of a dust 
component in the spectrum of the WC9 star can not be tested. 
It is however interesting to note that IRS13E2, the WN8 star 
analysed in Sect. I4.2.5I has a quite large K - L color. The entire 
cluster IRS 1 3E is also known to be at the top of a very prominent 
gas stream in L band dClenet et al.ll2004f) . Taken together, these 
arguments indicate that the stars of IRS13E may well all bathe 
in a continuum emission due to hot dust. We thus conclude that 
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Fig. 18. Best fit (red dot-dashed line) of the observed H+K band 
spectrum of IRS7SE (WC9, black solid line). 



their derived parameters must be regarded as only indicative (see 
also Sect. 16.51 ). 

With this restriction in mind, we show a tentative fit of the 
spectrum of the WC9 star IRS13E4 in Fig. [19] The main C lines 
are reasonably well reproduced by our model. However, the He i 
2.184 fj.m and Hen 2.189 fim lines, while having the right line 
ratio, are too strong. Reducing them would require a reduction 
of the mass loss rate. This is at the cost of a reduction of the 
C lines (which could be compensated by an increase in the C 
abundance), and of a modification of the shape of the lines. With 
lower M, lines get narrower and more centrally peaked because 
the density decreases. In conclusion, and given the above discus- 
sion, we argue that the near-IR spectrum of IRS13E4 is certainly 
contaminated by dust emission. Such a component could explain 
that the He i and He n lines are too weak (being diluted). Carbon 
lines are also certainly diluted, so that our C abundance estimate 
is likely a lower limit. We also probably over-estimate the total 
luminosity. 

4.4. OB supergiants 

Due to the rather limited signal to noise ratio of most of the spec- 
tra of the OB stars known in the central cluster (S/N of the order 
of a few tens on average), we have restricted ourselves to a gen- 
eral study of the properties of these stars. For tha t, we have used 
the ave rage spectra of 10 supergiants presented in Paumar d et alj 
(2006). Hence, we have derived average stellar and wind param- 
eters for this population of OB stars. 

Fig. l20l shows the best fit model. The main parameters used 
for this model are: T eff = 27500 K, logg = 3.25, log £ = 5.33, 

M = 3 x 10~ 7 M Q yr" 1 , v M = 1850 km s _1 , X(He) = 0.3 and v turb 
= 15 km s _1 . Additionally, a rotational velocity of ~ 100 km s _1 
could be estimated from the overall shape of all absorption lines 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




Mjitm) 

Fig. 19. Tentative fit (red dot-dashed line)of the K band spectra 
of the IRS13E4 (WC9, black solid line). See text for discussion 
of the discrepancies. 



0. The effective temperature is very difficult to constrain, just as 
in the Ofpe/WN9 stars studied previously. We can only put a 
constraint on the upper limit of T e ff from the absence of He n 
lines, especially Hen 2.189 ym. This upper limit is of around 
32000 K. Giving a lower limit is more challenging. We tried sev- 
eral values between 25000 and 32500 K. This range is thought 
to be appropriate since all stars contributing to our average spec- 
trum are late O / early B supergiants. For these t ypes of stars, 
T eff i s expected to be around 25000 - 30000 K (Martins et alj 
2005). Reasonable fits could be obtained for these different T e ff. 
As we have no diagnostic to better constrain T e ff, we finally 
adopted 27500 K as a typical value. For the luminosity, we also 
adopted lo g j- = 5 .3 since this is typical of late O supergiants 
(e.g. lMartins et al.ll2005l) . 

The mass loss rate is not strongly constrained either since 
most lines are seen in absorption and do not appear to be filled 
by wind emission. Here too, an upper limit on M of a few 10~ 6 
Mq yr~' can be given, above which Br y cannot be reproduced 
any more. 

Finally, the He abundance is tentatively constrained from the 
strength of He 1 2.1 12 pmi, He 1 2.150 ym, He 1 2.161 prni and Hei 
2.184 pun. The strength of these lines not only depends on the He 
content, but also on the microturbulent velocity. The slope of the 
velocity law (the (3 parameter) plays a role too. As we have no 
independent way of determining all these parameters, we tried 
different combinations with reasonable values (10 < v tur b < 20 
km s" 1 , 1 .0 < p < 2.0). In the end, we found that good fits could 
be achieved for He/H ratios in the range 0.2-0.35. This value 
will be discussed in Sect. [9] Finally, we stress that we do not 
reproduce the 2.115 fmx emission. As discussed in Sect. I4.1.T1 
this line is not clearly identified and we did not try to fit it. 



For that, the model spectra were convolved to take into account both 
the instrumental resolution and a rotational broadening represented by 
a simple Gaussian function. 
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Fig. 20. Best fit (red solid line) of the average spectrum of the 10 
supergiants (see lPaumard et alj|2006l) . 



5. Ionising radiation in the Galactic Center 

One of the crucial issues revealed by early studies of the Galactic 
Center was the apparent incompatibility between the ionisation 
of the gas and the ionising flux provi ded by the populati on of 
massive stars. This was highli ghted by Nai arro et al.l (1 19971) and 
more recently by iLutzl ( 1999|) in a detaile d analysis of ISO ob- 
servations (see also Thornle v et al.ll2000T) . In view of our new 
quantitative analysis, we argue that this issue is solved. Since 
this is an important result, we give a detailed explanation of the 
different aspects of the problem together with the proposed so- 
lutions. 



5.1. H and He i ionising photons 

The first part of the "ionisation problem" concerns the total num- 
ber of H and He i ionising photons produced by the GC massive 
stars (i.e. ionising photons short-ward of 912 A - Qu - and 504 
A- Quel-)- Nebular emission from the central H n region was 
used to constrain the various ionising fluxes. Radio m easure- 
ments of the free-free continuum bv lEkers et al.l (Il983l) showed 
that Qh = lO 50 - 4 * - 3 ^ 1 (see also Genz el etalJl!994h . Later, 
iKrabbe et al.l d!991l) derived Q Hel ~ lO^V 1 from a study of 
the nebular Hei 2.058 [im emission. These measurements could 
then be compared to the results of the qua ntitative modeling o f 
the stellar properties of the "He i" stars by Naia rro et al.l (Q997). 
Their conclusion was that the "He i" stars provided enough H 
ionising photons, but failed by several orders of magnitude to 
produce the extreme UV flux required to account for the nebular 
He ionisation. An underlying population of stars hotter than the 
"He i" stars but not detected was suspected to be responsible for 

this harder flux. 

Such a population has recently reported by Paumard et al. 
Hence, a re-estimate of the ionising power of the cen- 
tral cluster is required. Table [2] gives Qu and Q h e i for the stars 
analysed in the present paper. In addition to these, Paumard et al. 



(2006) identified 15 other evolved massive stars and 59 OB stars. 
To estimate the total ionisation flux delivered by this population 
we used the following approach: for th e WC9 stars not an alysed 
here, we adopted the ionising fluxes of Smith et al. (2002) using 
their model WC number 3; for the 2 Ofpe/WN9 stars IRS16NE 
and IRS16SW, we adopted the parameters of IRS16C; for the 
WN7 stars not studied, we adopted the values of AFNWNW, an- 
other WN7 star included in ou r sample; final l y, for t he O stars, 
we adopted the calibration of iMartins et a l. (2005) (see then- 
Table 4). Since all supergiants have a spectral type between 08.5 
and B2 but some of them suffer from a classification uncertainty 
of up to 2 spectral sub-types, we adopted the ionising flux of a 
09.51 star as typical of th e GC supergiants. As for O dwarfs, the 
values of Qu and Q Hel of lMartins et al.l d2005l) for a 09.5V star 
were chosen. For the B dwarfs, we simply adopted a value ten 
times smaller than for OV stars. This is a rough approximation 
which however has very small impact on the final result, B stars 
providing a negligible ionising flux. This leaves us with the fol- 
lowing numbers: Qu = 6.0 x 10 5 V and Q Hel = 2.3 x 10 4 V [ 
(or log Qu = 50.8 and log^Hei = 49.4). The contribution of 
the different classes of stars are gathered in Table [3] We thus 
conclude that not only the H ionising flux but also the He i ionis- 
ing flux required to reproduce the nebular emission can be pro- 
vided by the population of massive stars. In fact, they may even 
produce slightly more ionising photons, indicating that the H n 
region might be density bounde d. A final comment on the to- 
tal He i ionising flux is needed. Kra bbe et al.l (1 1 99 lb state that 
logfi H ei = 48.4 s"\ but also that Qu e i/Qu= 0.06. We find 
Qhsi/Qh- 0.04, in very good agreement. The difference in the 
absolute Qua values between our study and Kra bbe et al.l(ll99ll) 
is their lower reference Qu (log Qu - 49.6). 

At this point, a comment on the c ontribution of the "Hei" 
stars is necessary. Naia rro et al.l (Q997) showed that the 8 stars 
they analysed could account for most of the ionising radiation 
of the region. We see that with the current estimate, this con- 
clusion would not be valid. Why is that? The answer is rooted 
in 1) our lower bolometric luminosities and 2) the inclusion of 
line-blanketing in the atmosph ere models. This ing redient was 
not available at the time of the Naiarr oet al.l d 1997b study. Test 
models reveal that in such stars, line-blanketing effects lead to 
a large redistribution of the blocked UV flux to longer wave- 
lengths, mainly above 912 A. Consequently, Qu is reduced sig- 
nificantly. 



5.2. Ionising radiation and stellar population 

The second important issue concernin g the ionisi ng radiation in 
the Galactic Center was highlighted by Lutzl dl999l) . In his study 
of nebular fine structure mid-IR lines of metals observed with 
ISO, Lutz pointed out that stellar evolution appears to fail to ex- 
plain the GC massive stellar population. This claim was based on 
computations of population synthesis models for a single burst 
of star formation and a standard Salpeter IMF. After 7 Myrs, the 
age of the population thought to be appropriate at that time, the 
fraction of the total ionising luminosity provided by the part of 
the HR diagram where the "He i" stars are lying {logT stm < 4.5 
and L > 5.75 ) was of the order of 1 %. This was at odds with 
the results of Naia rro et alj d 1997b who argued that these stars 
could account for mo re than half the total ionising luminosity. 
This discrepancy lead lLuta (1999) to the conclusion that stellar 
evolution - indirectly tested here through synthesis population 
models - was not producing enough cool stars or equivalently 
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Type of star 


number ot stars 


Qh 
I s J 


fa total 


t^Hel 

I s J 


% t^Hel total 


OB V 


23 


3.1 x 10 4S 


0.5 


2.6 x 10 4b 


0.1 


OB I 


30 


3.0 x 10 50 


50.0 


6.7 x 10 48 


29.1 


Ofpe/WN9 


8 


2.5 x 10 49 


4.2 


3.4 x 10 48 


14.8 


WN 


9 


1.1 x 10 50 


18.3 


7.0 x 10 48 


30.4 


WN/C 


2 


2.0 x 10 49 


3.3 


1.1 x 10 48 


4.8 


wc 


13 


1.5 x 10 50 


25.0 


4.7 x 10 48 


20.4 


Total 




6.0 x lO 50 




2.3 x 10 49 





that the time spent by a massive star in the cool part of its track 
was much too short. 

This statement is no longer valid. Only six stars analysed 
here have logT stM < 4.5 and log j- > 5.75. And this is in the case 
we include IRS13E2 for which we only have an upper limit on 
log j-. To these four stars, we need to add the binary IRS16SW 
and the binary candidate IRS16NE. Both stars likely have prop- 
erties similar to IRS16C. In total, the H ionising flux of these 
stars represents about 9 % of the total Qu- This is an upper limit 
due to the possible overestimate of IRS13E2's luminosity. If we 
exclude IRS13E2, the remaining 7 stars contribute only 4 % of 
2^ tal - This is in excellent agreement with what is expected from 
a burst of star formation after 7 Myrs (w hich is within the ag e 
range now stated for the population, see iPaumard et al.l 12006). 
Once again, this is mainly due to the recent discovery of a hot 
population of OB and Wolf-Rayet stars responsible for the ma- 
jority of the ionising luminosity. Table [3] shows that the OB su- 
pergiants and Wolf-Rayet stars contribute more than 90 % of the 
ionising flux. The main conclusion is that standard stellar evolu- 
tion - used in population synthesis models - is able to account 
for the GC massive stars. 

5.3. Mid-IR nebular Ne lines 

The third problem with the ionisatio n of th e Gal actic Center 
region was also p ointed out by iLutzl ( 1 1999b (see also 
ThornlexetiD[2000). It concerned the low ionisation of the lo- 
cal gas as derived from the ratio of fine structure mid-IR lines of 
different ionisation states of Ne, namely Ne in 15.5 fim and Ne n 
12.8 jt/m. ISO observations revealed that [Ne m]/[Ne n] was 0.05. 
Using the SED predicted by the synthesis population model de- 
scribed in the previous section as an input of a photoionisatio n 
model performed with the code CLOUDY dFerland et alJI 19981) . 
Lutz revealed that after 7 Myrs, [Ne m]/[Ne n] was still as large 
as 1 to 2, or a factor 50 to 100 more than the observed value. 

To investigate this issue, we have performed photoionisa- 
tion models with CLOUDY (version C06.02). As an input SED, 
we have simply added all individual SEDs computed for the 
stars presently studied. For those stars whic h were not explicitly 
treated here, we have adopted the SEDs of Martins et al. (2005) 
or those of similar s tars present in the current sample. Adopting 
the same density as ILutzl d 1999b (3000 cm -3 ) and our total ion- 
ising flux implies an ionisation parameter log U = -0.60. With 
these values, we obtain a ratio [Ne m]/[Ne n] of 0.9- 1 .7 depend- 
ing on the geometry. This is still larger than the observed value 
(0.05). Note however that if we use the ionisation parameter as 



4 the ionisation parameter is defined by U = , 3p where r is the 
distance to the ionising source (chosen to be 0.5 pc in our case) and n is 
the density. 



ILutzl (1 19991) (lost/ = -l),we get [Nem]/[Ne n] ~ 0.5 -0.7. This 
is a factor 2-3 lower than the values of Lutz. 

How can we explain the still large values of [Nem]/[Neii]? 
One possibility is that we underestimate the density. There is ev- 
idence that values as large as 1 4 " 5 cm" 3 are required to produce 
[Oi], [Om] and [Fen] lines dGenzel et alJll984l Il985l) . Using 
such large values (and the corresponding ionisation parameter, 
log U = -1.2 and -1.6), CLOUDY models with our total stel- 
lar SED give [Nera]/[Neii] ~ 0.5 (for 10 4 cirT 3 ) and ~ 0.1 (for 
10 5 cm" 3 ). This is in better agreemen t with the observed ratio , 
although still a factor 2-10 too large. Shie lds & Ferlandl {1994) 
argued that the nebular spectrum of the Galactic Center could be 
reproduced only if several gas components with different den- 
sities were involved. In view of the present result, it may well 
be that the Ne ionisation requires a large density material (we 
can reproduce the observed [Ne m]/[Ne n] ratio for a density of 
3x 10 5 cm" 3 ). 

Another explanation to the large theoretical [Ne m]/[Ne n] 
ratio could be that we still overestimate the flux at 41eV, i.e. 
the Nen ionisation energy probed by the [Ne m]/[Ne n] ratio. 
This part of the spectrum is quite sensitive to blanketing ef- 
fects. A slight increase in metallicity could lead t o a reduced flux 
and c onsequently a lower [Ne m]/[Ne n] ratio (e.g. Morisse fet al.l 
12004 . 

One may also wonder which type of star contributes signif- 
icant flux at 41eV. Actually, it turns out that the total SED is 
completely dominated by a single star at this energy: the hot 
WN5/6 star IRS16SE2. To test the influence of this star on the 
ionisation of the GC gas, we removed its contribution to the to- 
tal SED and ran test CLOUDY models. Amazingly, for a den- 
sity of 3000 cm" 3 (and log U = -0.6), the [Ne m]/[Ne n] ratio 
drops to 0.008, less than the observed value! This shows that 
this ratio is extremely sensitive to the local radiation field. One 
can imagine that most of the ionised gas is not illuminated by 
the IRS16SE2 radiation due to shielding by local structures in 
molecular clouds. In that case, the remaining radiation field is 
soft enough to maintain a low [Ne m]/[Ne n] ratio. 

In conclusion, one could say that a revised nebular modelling 
taking into account both the spatial distribution of the gas and of 
the ionising sources is required to solve the Ne ionisation prob- 
lem. This is well beyond the scope of the present paper. 

6. Stellar evolution in the Galactic Center 

In this Section we discuss, in view of the results of our quanti- 
tative analysis, the evolution of massive stars beyond the main 
sequence. It is generally accepted that stars in the mass range 
25-60 M Q evolve from O stars to WN H-poor stars through a 
LBV and/or red supergiant phase before becoming WC stars (for 
M > 40 M ). In the following, we refine this scenario in the 
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particular case of the GC, establishing a plausible evolutionary 
sequence between Ofpe/WN9, WN8 and WN/C stars (Sect.loTTI 
16.2116.31 ). We also quantitatively compare the position of the GC 
Wolf-Rayet stars in the HR diagram to evolutionary tracks in- 
cluding rotation (Sect. [6~4b . Finally we discuss the properties of 
the IRS13E cluster stars dOt . 



6.1. Ofpe/WN9 stars 

In the present study, we have ana lysed five stars classified 
Ofpe/WN9 by iPaumard et all 12006). Amon g these five stars, 
three have been studied by iNajarro et all dl997l) : IRS16C, 
IRS16NW and AF. Compared to the Najarro et al. analysis, we 
find a similar range of luminosities (although for IRS16C our 
luminosity is lower) and the same terminal velocities (within the 
uncertainty). The effective temperatures are higher in our study, 
partly due to the inclusion of line-blanketing in our models as 
already discussed. Nevertheless, T e ff remains poorly constrained 
so that th e range of acceptabl e values overlap with the temper- 
atures of Naja rro et al.l (jl997). As a consequence of the hotter 
T e ff, our radii are smaller. But the main differences concern 1) 
the mass loss rates and 2) the He content. Both parameters are 
linked to some extent: when fitting Hei 2.058 fim, Hei 2.112 
[im and Br r , adopting a larger He/H content will require a larger 
M in order to reproduce the level of Br y emission. Of course, 
in that case Hei 2.058 fim and Hei 2.112 /urn get stronger too. 
But their absolute strength is also controlled by the He ionisa- 
tion which in turn depends on the effective temperature and the 
line-blanketing effect. Since we used more realistic atmosphere 
models as well as better spectra (higher S/N ratio and spectral 
resolution, good correction from nebular emission), we argue 
that our d erived parameters r epresent an improvement over the 
result of Naja rro et al.l d 1997b . In practice, we find values of M 
3 to 10 times lower than Najarro et al., and much lower He con- 
tents (He/H ~ 0.2-0.5 compared to 1.3-3.0). 

These revised parameters are important for assessing the 
evolutionary status of the Ofpe/WN9 stars (see next Section). 
They are also very interesting since they bring the GC 
Ofpe/WN 9 stars closer to other Galactic and LMC stars of 
this type. ICrowther & Smith! (Il997l) analysed a sample of LMC 
Ofpe/WN9 stars and found that their He content was much 
smaller than in the GC stars, which was tentatively attributed 
to possible metallicity effects. Our new values are in better 
agreement with the Crowther et al. measurement, showing that 
the LMC and GC Ofpe / WN9 stars are chemically similar. In 
contrast, iPasquali etall d 19971) found He/H ~ 0.5 (by num- 
ber) for a sample of LMC Ofpe/WN9 stars partly overlapping 
with the Crowther et al. sample. These He contents are only 
marginall y larger than ours, and are certainly lower than the 
values o f Naja rro et alj dl997l) . Concerning the mass loss rates, 
we find that on avera g e M is sys tematically smaller t han in 
the ICrowther & Smith! d 19971) and IPasquali et all d 19971) stud- 
ies. This difference is surprising since the lower metallicity of 
the LMC should lead to lower mass loss rates. Combined with 
the slightly larger He content, this may be an indication of a 
more advanced evolution. Another indicator points to the same 
conclusion: the terminal velocities of the GC Ofpe/WN9 stars 
are usually larger than for stars of the s ame spectral type (see 
Pasq uali et aD 19971: iBresolin et al.ll2002b . In terestingly enough, 
the on ly Galactic star of this type studied by ICrowther & Smith! 
1997) had a much larger Voo than the LMC stars. ISchaererl 



1996) explained such a trend by a more advanced evolution in 



the Ofpe/WN9 phase for stars in the LMC, with the consequence 



of greater proximity to the Eddington limit and implying a lower 
terminal velocity. 

In conclusion, in view of our study, the class of Ofpe/WN9 
seems to be more homogeneous than previously believed. The 
difference between GC and LMC Ofpe/WN9 stars is likely due 
to a different state of chemical evolution: massive stars become 
Ofpe/WN9 stars later in a low Z environment such as the LMC, 
mainly due to lower mass loss rates. 

6.2. An evolutionary link between Ofpe/WN9 and WN8 stars 

On the basis of their K band morphology, all of the initial "He i" 
stars have been class ified as either Ofpe/WN9 or WN8 stars by 
Paumard et al. (2006). These two spectral types are indeed char- 
acteristic of relatively cool evolved massive stars. In the K band, 
they are defined by strong He i and Br y emission lines. The rela- 
tive intensity of these different lines is similar in both spectral 
types: Hei 2.058 /urn is usually stronger than Br y , while Hei 
2.112 /urn is weaker or of equal strength. The main differences 
are 1) the shape of the lines (Ofpe/WN9 stars show P-Cygni pro- 
files in the He i lines, while WN8 stars have strong pure emission 
lines), 2) their width (WN8 stars have broader lines), 3) their ab- 
solute strength (stronger lines in WN8 stars) and 4) the presence 
of weak Hen lines in WN8 stars. Inspection of Table [2] reveals 
that quantitatively, these morphological differences are due to 
larger mass loss rates (and to a lesser extent wind terminal ve- 
locities) as well as higher effective temperatures for WN8 stars. 
Indeed, T eff ranges from 30000 to 41000 K for WN8 stars, while 
they are lower than 30000 K for Ofpe/WN9 stars. Mass loss rates 
are ~ 2-4 times smaller in Ofpe/WN9 stars. Luminosities are 
also slightly lower in WN8 stars. This comparison indicates that 
Ofpe/WN9 stars and WN8 stars may be physically related and 
may well represent consecutive phases of a single evolutionary 
sequence. 

Fig. |2T] displays the H content as a function of luminosity 
in evolutionary models (solid lines) and shows the position of 
the Ofpe/WN9 and WN8 stars. It is clear that both types of 
stars gather in different parts of the diagram: Ofpe/WN9 stars 
still show a significant amount of hydrogen in their atmospheres 
whereas WN8 stars are mainly H free. This can be interpreted 
as an evolutionary sequence where Ofpe/WN9 stars evolve into 
WN8 stars. Such a scenario would be consistent with the prop- 
erties reported above. Ofpe/WN9 stars could well be on the cool 
part of an evolutionary track. This track will then loop back to 
the hot part of the HR diagram. As a star follows this track, it 
evolves chemically, gets hotter and strengthen its wind on its 
way to the WR phase. The hotter T e ff and larger M explains the 
appearance of Hen lines and the stronger emission lines, and the 
lower H content reveals the chemical evolution of the star. In this 
scenario, the AF star could be in an intermediate state. Its spec- 
tral morphology is similar to WN8 stars except that it does not 
show Hen lines, indicating a relatively cool T e ff (confirmed by 
the quantitative analysis, see Table |2j. Its mass loss rate is also 
more typical of WN8 stars. In Fig.|2T| the AF star lies in between 
the groups of Ofpe/WN9 and WN8 stars, with a H mass fraction 
< 0.3. Hence, it is also more evolved than Ofpe/WN9 stars, but 
less than WN8 stars, and nicely fits in the evolutionary scenario 
we suggest. 

Another argument in favour of this scenario is the variabil- 
ity of both types of objects. On the one hand, the Ofpe/WN9 
stars in the Galactic Center have been claim ed to be LBV candi- 
dates, or even LB Vs in a quiescent phase (Pau mard et alj |2004: 
Trip pe et alj2006l) . This is based on spectral similarities between 
these stars and objects elsewhere in the Galaxy and LMC known 
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Fig. 21. Hydrogen mass fraction as a function of luminosity. 
Solid lin es are the Geneva evoluti onary tracks with rotation and 
Z = Z dMevnet & MaedeJl2005b . The ZAMS masses for each 
track are indicated. Stars show the position of the Ofpe/WN9 
stars (IRS16C, IRS16NW, IRS33SE, IRS34W, AF) analysed 
here, while the triangles are the WN8 stars of our sample. The 
AF star lies at log ^ = 5.3 and X(H)=0.3. See Sect.|6~4]for dis- 
cussion. 



to be related to LBV stars. Besides, one of them - IRS34W 
- was shown to be photometrica lly variable on timescales of 
months-years. Trip pe et"aT] J2006) interpreted that as a sign of 
obscuration by dust produced in material ejected in an LB V-type 
outburst of the star. On the other hand, WN8 stars are known 
to be the c lass of Wolf-Rayet stars experiencing the str ongest 
variability dAntokhin et all 1 19951: iMarchenko et alJll998l) . This 
high degree of variability may be related to the LBV phe- 
nomenon. A link between WN8 stars and LBVs is also favoured 
by the presence of LB V-like nebulosities around most of them 
dCrowther et alj|1995bl) . 

We thus argue that the GC Ofpe/WN9 stars are precur- 
sors of WN8 stars, and are in a state closely related to the 
LBV phase. This pictur e is fully consistent with the scenario of 
ICrowther et alj d!995bl) . further extended bv lCrowther & S mith 
1997b : a 25-60 Mo star evolves into a WN9-11 star (similar 



to Ofpe/WN9, see lCrowther et al]|l995al) before experiencing a 
LBV phase and becoming a WN8 Wo lf-Rayet star. T h e prop- 
ertie s of the WN8 star s analysed by ICrowther et alj dj995b) 
and iHerald et al.l d200ll) are very similar to those of our sam- 
ple WN8 stars. The only difference is a larger s pread in H con- 
tent: wh ile we find all WN8 to be almost H free. ICrowther et all 
dl995bl) have both H fre e stars and stars still showing a signifi- 
cant amount of hvdrogen. lHerald et al.l d2001l) also found X(H) ~ 
20% by mass. However, this H content remains much lower than 
in the GC Ofpe/WN9 stars. Hence, the suggested evolutionary 
scenario remains valid. 



6.3. WN/C stars 

WN/ C stars are Wolf-Rayet stars showing both strong C and N 
lines dMassev & Grovd[T989HWillis & Stickland|[T99a) . Carbon 
is produced by the triple a reaction, while Nitrogen results 
mainly from CNO processing. It is widely accepted that WN/C 
stars are core He burning stars with a CNO-enriched envelope 
in which mixing processes have created a layer with both H and 
He burning products. W hen mass loss reveals this layer, the star 
turns into a WN/C star (lLangeril 19911 iMeynet & Maederll2003l 
120051) . Quantitatively, such WN/C stars have C/N ratios of the 
order 1 . 

For a long time, the fraction of WR stars in the WN/C state 
observed in the Galaxy has been difficult to explain with evo- 
lutionary models. Abundance profiles in such models were usu- 
ally very steep in the transition region between H and He burn- 
ing products, so that the region where both type of products 
were present was extremely thin. Consequently, it was quickly 
removed by the stellar wind, resulting in a very short lifetime 
of the WN/C phase. Consequently, the number of WN/C stars 
predicted by such models was much lower than the observed 
value. Significant improvements have been made in the last years 
mainl y due to the inclusion of mi xing processes triggered by ro- 
tation. iMevnet & Maederl d2003l) have shown that rotation cre- 
ated shallower abundance gradients in stellar interiors, increas- 
ing the size of the mixed H and He burning products. As a conse- 
quence, the lifetime of the WN/C phase is lengthened, resulting 
in a total number of WN/C stars in better agreement with obser- 
vations. 

What about the GC Wolf-Rayet population? Our analysis of 
IRS15SW and IRS7SW has revealed that they were significantly 
enriched in Carbon compared to other WN stars, leading to a 
classification as WN8/WC9 stars. Fig. l22l shows the position of 
our sample stars in a log C/N - logC/He diagram. It turns out 
that IRS15SW and IRS7SW lie in between WN stars (stars with 
low C content) and the C-rich WC9 stars. Their position is also 
in excellent agreement wit h WR8 and WR145, two WN/C stars 
studied by ICrowther et al.l (Il995dl) . and with the prediction of 
evolutionary models. Hence we have a quantitative confirmation 
that IRS 15SW and I RS7SW are core He burning objects on then- 
way to a WC phase. IMevnet & Maederl d2003l) argue that only 
stars with initial masses in the range 30-60 M go through the 
WN/C phase: more massive stars have very strong winds that 
quickly remove the CNO enriched envelope; lower mass stars 
have a too H-rich envelope and He burning products are too 
diluted. Estimates of the present-day masses of IRS15SW and 
IRS7S W from the mass luminosity relation of Heger & Langer 
(1 19961) for H free WR stars gives 10.3 and 20.0 M Q respectively, 
implying that these stars have lost 30 to 80 % of their mass 
through stellar winds. 

If we compare the properties of IRS15SW and IRS7SW to 
WR8 and WR145 (see ICrowther et al.|[l995cJ) . the luminosities 
are all similar, in the range 10 51 ~ 55 L Q , which in turn implies 
that the present-day mass of these four objects are also very 
close (10 - 20 M ). The clumping corrected mass loss rates 
are also similar, with logM/ ^[f ~ -4.3. However, the ter- 
minal velocities of the GC stars are smaller than for WR8 and 
WR145 (700-800 km s _1 as opposed to 1390-1590 km s" 1 ) as 
well as the stellar temperatures (32-37 kK vs 41-48 kK). These 
differences reflect the different spectral types of the two samples: 
while the GC WN/C stars are late type WR star s (WN8/WC9) 
WR8 and WR145 are earlier (WN6/WC4, see ICrowther et ail 
1995d). Earlier WN and WC type stars have higher effective 
temperatures and larger terminal velocities (see Fig. |25l . This 
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can explain the observed trend: both types of stars (early and late 
WN/C stars) have the same luminosity but early WN/C stars are 
hotter, which implies that their radius is smaller. Consequently, 
their escape velocity, scaling as {M/R) 5 , is larger (the present 
mass being approximately the same). Since the terminal velocity 
is directly proportional to the escape velocity, one naturally finds 
larger in early WN/C stars. 

We argue that the GC WN/C stars are most likely the de- 
scendents of WN8 stars. Their spectral morphology is extremely 
similar, except for the presence of Civ and Cm lines in the 
WN/C stars. The quantitative analysis confirms their "twin" 
character: the ranges of values for T e ff luminosities, mass loss 
rates, terminal velocities, He and N content are the same for 
both WN8 and WN/C stars (see Table |2). The only possible 
exceptions are IRS13E2 (but see discussion in Sect. 16. 5\ and 
IRS9W. Fig.|22]shows that the WN/C stars are more chemically 
evolved than WN stars. Hence, we conclude that most of the 
GC WN8 stars will go through a WN/C phase as two of them 
are currently doing. 

Summarizing the last three Sections, we argue that in the 
Galactic Center, stars with initial masses in the range 30 - 60 
M follow the evolutionary sequence 

(Ofpe/WN9 ^ LBV) -> WN8 -> WN/C 

on their way to the supernova explosion. On the main sequence, 
they probably appear as mid/early O stars. After the WN/C 
phase, they most likely become late WC stars. Indeed the ef- 
fective temperatures of WN/C stars are slightly lower than the 
WC9 stars of our sample. Hence, we can expect them to enter the 
WC sequence from the low ionisation (or equivalently the low 
T eff ) side and appear as WC9 stars. The (Ofpe/WN9 ^ LBV) 
sequence in the suggested scenario indicates that Ofpe/WN9 and 
LBV stars are closely inter-related and most likely represent dif- 
ferent states of a single evolutionary phase (see discussion in 
iTrippe et alj|2006h . 

6.4. Are stellar evolutionary tracks too luminous? 

The direct comparison between derived parameters and evolu- 
tionary tracks by means of a classical HR diagram (logL - 
logT e ff) is not straightforward for Wolf-Rayet stars. In evolu- 
tionary models, the "effective" temperature (noted T evo i in the 
present discussion) is defined at the outer boundary by Eq. 13.21 
However, there is no atmosphere in such models and strictly 
speaking, the radius in this case (noted R eV oi) is not the one 
for which the optical depth is equal to 2/3 (exact definition of 
T e ff). It corresponds more to a hydrostatic radius. In Wolf-Rayet 
stars, such a radius - for which the atmosphere is quasi-static 
- is at an optical depth much larger than 2/3. Said differently, 
R2/3 for which T e ff can be defined is much larger than R evo i- 
Consequently, T e ft is lower than T evo i- To overcome this problem, 
one usually defines a radius R* (and the corresponding temper- 
ature T*, see Eq. \3.2\ at an optical depth equal to 20. Such a 
radius is more comparable to R evo i so that T* and T evo i can be 
directly compared. 

In Fig. [23] we show the HR diagram built using T* for our 
program stars (See Table [2]). They are shown by the filled sym- 
bols (see next Section for a discussion of the two stars repre- 
sented by open symbols). The evo lutionary tracks including ro- 
tation o flMevnet & Maederl (12005) for solar metallicity (left) and 
twice solar metallicity (right) are overplotted. In the diagram at 
Z = Z Q , we clearly see that the Ofpe/WN9 stars populate the 




log C/He (by number) 



Fig. 22. logC/N as a function of log C/He (by number) 
from evolutionary tracks with rotation (dot-dashed line) from 
iMevnet & Maederl d2005l) . Filled symbols are the WR stars anal- 
ysed in the present paper (triangles: WN stars; diamonds: WN/C 
stars; circles: WC9 stars ). Open diamonds are t he WN/C stars 
WR8 and WR145 from (Crowther et ail dl995dl) . IRS7SW and 
IRS15SW, the two WN/C stars of our sample, have C/N and 
C/He ratios very similar to WR8 and WR145. The GC WN stars 
all lie at log C/He < -3. For these stars, we have only upper 
limits on the C abundance, and consequently on log C/N and 
log C/He. The WC9 stars analysed here have large C/He ratios 
and only lower limits on C/N since no N lines are present in 
the spectra and only an upper limit on the N abundances can be 
estimated. 



coolest region, while the more evolved Wolf-Rayet stars are hot- 
ter and, on average, less luminous. Qualitatively, this trend is 
similar to the 40 M track: coming back from its redward ex- 
tension, this track goes to lower temperatures. However, this 
track, as well as the ones for other masses, always remains at 
log L/L G > 5.4, in contrast to most of the Wolf-Rayet stars of 
our sample which have 5.1 < log L/L e < 5.5. Hence, we see 
that quantitatively the solar metallicity evolutionary tracks do 
not seem to explain the evolution of the GC Wolf-Rayet stars. 
Note that this is independent of any definition of the tempera- 
ture, since the problem is due to luminosities. 

If we now focus on the HR diagram where the Z = 2 x Z Q 
metallicity tracks are plotted, we see that lower luminosities are 
reached. This is due to the stronger mass loss rates in these 
tracks which "peel off" the star more quickly, reducing its ra- 
dius and consequently its luminosity. Does that mean that the 
GC stars have a supersolar metallicity? We will see in Sect. 
Q that the question is still largely debated. The position of the 
GC Wolf-Rayet stars in the HR diagram may be an indica- 
tion of a super-solar environment. But this can also be due to 
inadequate mass loss rates in evolutionary calculation at solar 
metallicity. If the amount of mass lost during the Wolf-Rayet 
phase is underestimated in such calculations, the tracks will be 
too luminous. Large M can be produced by a fast rotation (see 
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iMaeder & Mevne3 l2QOO). but the analysis of OB stars average 
spectrum do not point to any particularly high V sin/ (see Sect. 
[9}. As long as the GC metallicity remains poorly constrained, a 
quantitative test of evolutionary tracks and their input parame- 
ters, such as mass loss rates, is thus not feasible. We will see in 
Sect. |7]that more work is needed to narrow the range of accept- 
able values for Z in the Galactic Center. 



6.5. The stars oflRS13E 

In the previous Section, we have deliberately excluded two stars 
from the discussion: the ones represented by open symbols in 
Fig. [23] These stars are the two Wolf-Rayet components of the 
cluster IRS13E, namely IRS13E2 and IRS13E4. The reason is 
that they show peculiar properties compared to the bulk of the 
GC Wolf-Rayet stars. Take first IRS13E2. It is a WN8 star (see 
its spectrum in Fig.fTTTi with parameters in marginal agreement 
with the other WN8 stars of our sample, with the notable ex- 
ception of its luminosity (log L/L Q = 6.1) which is much larger 
than the range of value s for WN8 stars (5.1-5.5, see Tableland 
ICrowther etaITl 995b). The same statement is true for IRS13E4. 
It is a WC9 star but exceptionally bright: whil e most stars of 
this ty pe have log L/L e = 5.0 (see an example in Crowther et al. 
2006), IRS13E4 is nearly 10 times more luminous. In addition, 
it shows a large, although not unprecedented, terminal velocity 
(see Fig.l25ll. 

How can we interpret these results? First, the most obvi- 
ous explanatio n is binarity. IRS 13E was initially thought to be 
a single star (Rrabbe et al. 1991) before being resolved into sev- 
eral components dMaillard et al .1120041) . The cluster is extremely 
compact and it is quite possible that the bright components might 
actually be multiple stars. But even with better spatially re- 
duced images, it may not be sufficient to distinguish between 
the components of a massive binary. Spectrophotometric moni- 
toring should then be used to test the binary hypothesis. So far, 
IRS13E2 and IRS13E4 (Fig. \U\ and [19]) do not show any sign 
of double line spectrum. In addition, photometric observations 
over the recent years (but with a very coarse sampling) show a 
flat light-curve for both stars (Trippe et al., in prep.). Better mon- 
itoring is required, but so far no clear evidence for binarity exists 
for IRS 1 3E2 and IRS 1 3En. 

But the most important factor was mentioned in Sect. 14.31 
dust contamination. There is evidence of a large dust content 
in IRS13E which may affect the infrared spectra of its stel- 
lar components. As a consequence, the luminosities are most 
likely upper limits. However, the luminosity difference between 
IRS13E4 and IRS7W is 0.7 dex, or a factor of 5. Assuming 
that both stars have the same luminosity and same spectral en- 
ergy distribution, this means that the total K band flux is 5 
times lar ger than the stel l ar+win d flux. This is exactly what is 
found by ICrowther et alj d2006l) in their analysis of the WC9 
star HD 164270: dust contributes 80% of the continuum. Note 
that their estimate is at 3 /im and not in the K - band, where dust 
contamination is smaller. But ICrowther et al] ([2006 ) also claim 
that HD 164270 has a weak dust shell compared to other dusty 
WC stars. Hence, for a typical WC star a 80 % dust contribution 
to the K band continuum is not unrealistic. 

Although further high spatial resolution photometric data in 
the near/mid IR range are needed to establish the true SED of 
each component, it is likely that dust explains most of the pecu- 
liar properties of the IRS13E Wolf-Rayet stars. 



7. Metallicity in the central parsec 

We first present a determination of the global metallicity in the 
GC by means of WN stars, and we then discuss its effect on He 1 
2.058 ymi. 

7. 1 . Metallicity from N content of WN stars 

Naia rro et al.l d2004h presented a new method for deriving metal- 
licity in evolved massive stars. Their idea relies on the fact that 
the surface Nitrogen content of massive stars reaches a maxi- 
mum during evolution due to production of N through the CNO 
cycle. Observationally, this maximum is obtained in WN stars. 
Its exact value is independent of the initial mass and of the wind 
properties. It changes with the initial metal content since the 
amount of Nitrogen synthesized is directly linked to the initial 
CNO content. Hence, the comparison of derived N abundances 
in WN stars to evolutionary tracks is an indirect tracer of the 
initial metallicity. 

In Table|2]we list our derived Nitrogen mass fraction. Only 
for WN stars such an estimate could be performed due to the 
presence of N m lines in the K band spectrum, especially N m 
2.247, 2.251 nm. For most stars, we have X(N) ~ 0.0135, with 
the exception of IRS15SW and AFNW which show a larger N 
content (up to 0.0326). However, for the latter star the uncer- 
tainty on the N abundance determination is quite large due to 
the low resolution and low S/N ratio of our spectrum. A direct 
comparison of this range of values to predictions of surface en- 
richm ent in evolutionary models (see e.g. Fig. 4 of Naiarro et al. 
12004 indicates an initial metallicity between solar and twice so- 
lar. If IRS15SW and AFNW are not taken into account, then a 
solar metallicity is favoured. Strictly speaking, we derive a lower 
limit since one cannot be sure that all WN stars have reached the 
phase in which their surface N abundance is maximum. Indeed, 
although this phase is long for initially very massive stars, it is 
quite short for stars with M ~ 20-50 M Q which is more appro- 
priate for our WN8 stars. In conclusion, we estimate the stellar 
GC metallicity to be at least solar. Stronger constraints cannot 
be derived from the present set of data/models. 

Stellar studies usuall y indicate a solar metallicity for the 
GC. iNaiarro et all d2004l) found Z = Z© for WN stars in the 
Arches cluster. T his i s also comparable to th e determinations of 
ICarr et al] d2000t) and lRanurez et al] ([2000) who also derived a 
solar [Fe/H] from the study of Iron lines in red supergiants in the 
central 2.5 pc. 

Abun dances have also been derived from interstellar gas 
studies. Shields & Ferlandl 0994) derived solar Ne content and 
twice solar Ar and N based on photoionisation models aimed 
at reproducing nebular fine structure lines in the mid infrared. 
However, improvements in the knowledge of the Ar collisional 
strengths led to a revision of th e Ar abundance d own to a nearly 
solar value (see discussion in ICarr et aL 2000). This is to be 
contra sted by the recent analysis of ISO data by iGiveon et al] 
d2002l) who found Ne ~ 1.4 N e H and Ar ~ 2.5 Ar Q . Similarly, 
Martin-Herna ndez et al] (2003) revised the Galactic metallic- 
ity gradients. Extrapolating their results to the Galactic Center, 
one should expect n(Ne/H) ~ n(Ar/H) ~ 1.5 - 2.5Z and 
n(N/H) ~ 5Z . 

Finally, o bservations o f K sh ell Fe lines by Chandra have 
been used by Mae da et~ai1 d2002) to derive an Iron constant as 
large as 4 times solar for the H 11 region SgrA East. Such a large 
metallicity was recently c onfirmed on a larger scale by X-ray 
observations with Suzaku (Kovama 20Q6|). 
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Fig. 23. HR diagram of WR stars in the GC. Left: solar metallicity tracks with rotation (V sin/ = 300 km s on the main sequence). 
Right: twice solar metallicity tracks including rotation. Star symbols are Ofpe/WN9 stars, triangles WN8 stars, pentagons WN5-7 
stars and circles WC9 stars. Different line thickness indicates different evolutionary tracks (the thicker the line, the higher the mass). 
ZAMS masses are marked for each track. The Humphreys-Davidson limit is also shown in the right upper part of each diagram. 



Clearly, the measurement of the metallicity in the Galactic 
Center needs more investigation. We have seen in Sect. 14.1.31 and 
I4.2.1l that individual elemental abundances could be super-solar 
(Mg, Si). But the uncertainties associated to these estimates did 
not allow a reliable statement. The use of high signal to noise, 
high spectral resolution data is mandatory in order to resolve 
weak metallic lines of Ofpe/WN9 stars, red supergiants as well 
as AGB stars. The distinction between a element and Iron abun- 
dances is also crucial, a elements are essentially produced in 
massive stars while Iron is mainly synthesized in low mass stars 
during the type la supernova phase. The ratio of a elements to 
Iron abundances is thus an indirect tracer of the slope of the IMF. 



7.2. Metallicity effect on He i 2.058 ymi 



We have seen in Sect. 14.2.31 that Hei 2.058 /urn is extremely sen- 
sitive to the metal content. Indeed, the higher the metallicity, 
the stronger the line emission so that Z = 2 x Z Q was actu- 
all y preferred for IRS9W . Su ch a behavior was previously n oted 
bv lCrowther et al.l d!998h and lBohannan & Crowtheri d 1999b and 
was attributed to the sensitivity of He i 2.058 fjm to the EUV ra- 
diation field, which in turn depends critically on the heavy-metal 
content. Fig.|24]shows the He ionisation and temperature struc- 
ture of the models for IRS9W with Z = Z and Z = 2 x Z Q . This 
figure helps us understand in detail what happens when the metal 
content is increased. Due to the larger blocking of radiation by 
metal opacities, the EUV flux is reduced which reduces the ioni- 
sation in the outer part of the atmosphere (see lower panel of Fig. 
|24"1 >. In addition, the escape of radiation through metallic lines is 
favored so that the efficiency of line cooling is enhanced, leading 
to a reduction of the temperature for log TR OSSe iand ^ - 1 .0 (upper 
panel of Fig. l24l i. Both effects (line cooling and line blocking) 
increase the Hei content in the outer atmosphere. In Fig. [24] we 
also indicate the line formation region of He 1 2.058 ymi and He n 



2.189 pmi. We see that He n 2. 1 89 fan emerges from a zone where 
the atmospheric structure is barely modified, so that its strength 
is similar in both models (see Fig. |9). On the other hand, He i 
2.058 fim is produced in the outer atmosphere where the amount 
of He i is larger in the Z = 2 x Z model. Consequently, He i 
2.058 fxm is stronger in the high metallicity model. Note that this 
explan ation is fully consistent with the approach of Naiarro et al. 
(119941) : when metallicity is increased, the He i content increases 
which leads to a decrease of the escape probability of the He i 
584 A line controlling the upper level of He 1 2.058 ym - He n be- 
ing the dominant ionisation s tate. This implies a stronger emis- 
sion at 2.058 //m(see Eq. 2 of lNaiarro etalll 19941) . 

Does that mean that a twice solar metallicity should be 
adopted for IRS9W? This would be an over-interpretation of the 
results. Indeed, what the previous exercise reveals is that Hei 
2.058 pmi is extremely sensitive to the EUV radiation. We cannot 
be absolutely sure that we correctly predict this part of the spec- 
trum. Indeed, our models are limited in the sense that we cannot 
include all lines from all elements. Hei 2.058 fim may also suf- 
fer from inaccuracies in metal atomic data. Hence, it may well 
be that increasing the metal content artificially compensates for 
the lack of metallic lines / elements. Besides, we have discussed 
only one example for which a larger metal content improves the 
fit of He i 2.058 /urn. But there are other stars (e.g. IRS7E2) for 
which having a higher Z does not help, Hei 2.058 /urn remain- 
ing too weak. Hence, we conclude that the metallicity cannot be 
accurately derived from its effect on He i 2.058 ymi. 



8. Wind properties of Wolf-Rayet stars 

We have derived the wind properties of a homogeneous and rea- 
sonably large sample of Wolf-Rayet stars, mainly composed of 
WN and Ofpe/WN9 stars. In order to see if the GC stellar winds 
are similar to other Galactic WR stars, we inspect the behav- 
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1°2 T Rosseland 

Fig. 24. Temperature (upper panel) and He ionisation (lower 
panel) structure of the best fit model for IRS9W. The solid line 
is the model with solar metallicity, while the dashed line is for 
Z = 2 x Z Q . In the upper panel, the line formation regions of He i 
2.058 fan and He n 2. 189 pmi are indicated by horizontal lines. 
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Fig. 25. Terminal velocity of GC WR stars (filled symbols) and 
other Galactic stars (open symbols) as a function of spectral 
ty pe. The data fo r galac tic stars are taken from the compilation 
of van der Huchtl (|200l[). Upper panel: WN stars; Lower panel: 
WC stars. 



ior of the terminal velocities and mass loss rates (corrected for 
clumping). 

Termi nal velocities fo r other Galactic stars are taken from the 
catalog of van der Huchtl d200lh and are compared to our sam- 
ple in Fig. [25] We see that there is a remarkable agreement. We 
have seen in Sect. 16. II that Voo is usually larger in the GC than in 
the LMC for Ofpe/WN9 stars, but this was at least qualitatively 
explained by a different metallicity/evolutionary state. We then 
conclude that the terminal velocities of the GC Wolf-Rayet stars 
are identical to other Galactic stars. 

The mass loss rates of WN and Ofpe/WN9 stars are shown 
in Fig. |26l as a function o f luminosity q We also p l ot the results 



of Crowth eTetaTI dl995d) . ICrowfher et al1jl995bl).lMorris et alj 
(120001) . (all in blue open symbols) and lHamann et al.l (2006) 
(black filled symbols). The GC mass loss rates are in good agree- 
ment with the Crowther et al./Morris et al. results for late WN 
stars (triangles). For early WN stars, the samples are too small 
to draw any conclusion. On the contrar y, there is a large dis crep- 
ancy between our results and those of Hamann et al. (2006|), at 
least for late WN stars. This difference is not so much due to 
the mass loss rates themselves as it is to luminosities. Indeed, 
both samples have similar values of M@J Jf, but the lumi- 
nosities in the GC are lower. Since most of the stars of the 
Hamann et al. sample have poorly constrained distances (only 
a few are in clusters), while the distance to the GC is accu- 
rately known, we argue that the main reason for the observed 
difference is an overestimate of the luminosities in the Hamann 
et al. sample. An additional explanation to the GC/Hamann et 
al. WNL luminosity difference could be that their sample con- 
tains several H-rich WN stars. These stars are thought to be 
very massive stars with a Wolf-Rayet appearance but possibly 



still core H burnin g objects on/close to the main sequence (e.g. 
Moff at et all 120061) . Such stars are expected to be very lumi- 
nous. How ever, they usually have spe ctral subtypes earlier than 
WN7 (e.g. lCrowfher & DessartlfT998l) . A significant fraction of 
the very luminous stars in the Hamann et al. sample (represented 
by filled triangles around log ^ ~ 6.0-6.2 in Fig. [26l> do not fall 
into this category (they are later WN stars). Their luminosity is 
thus most likely overestimated due to poorly known distance. 

In Fig. [26] we also plot the mass loss rates of the GC 
Ofpe/WN9 stars together with Galactic O supergiants. We see 
that in terms of M, the Ofpe/WN9 stars are closer to O stars than 
to WN stars. This is another indication that Ofpe/WN9 stars are 
precursors of WN8 stars as we have argued in Sect. 16.21 



9. He abundance in OB stars 



5 our sample of WC stars is too small for a comparison to other 
Galactic stars to be relevant 



We have seen in Sect. 14.41 that the stellar and wind parameters 
of the average OB stars in the central parsec were not strongly 
constrained. Most of them have only upper limits (TW, M). For 
the helium content, we could estimate X(He) ^ 0.2 - 0.35. This 
value is larger than the standard solar abundance (0.1). 

Since OB supergiants are massive stars evolving away from 
the main sequence, chemical enrichment is the most likely ex- 
planation. The question is whether or not the amount of He ob- 
served is consistent with such a mechanism. In that respect, the 
inclusion of rotation in evolutionary models is a key ingredient 
since it triggers additional mixing which mod ifies the He surface 
abund ance X(He). In the recent models of iMevnet & Maederl 
(2005), X(He) evolves slowly from the initial value during all the 
"O" phase, and then dramatically increases when the star enters 
the Wolf-Rayet phase. At the end of the O phase, a He content as 
large as 0.3-0.4 can be obtained. This is compatible with our de- 
rived value. We have also seen in Sect. 14. 41 that the average pro- 
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Fig. 26. Mass loss rate corrected for clumping (M/ yjf) as 
a function of luminosity for GC WN stars (red filled sym- 
bols) compared to other Galactic Wolf-Rayet stars (blue and 
black symbols) and O supergiants (magenta star symbols). 
Compa ri son data from ICrowther et al] 



[(ll995bl).ICrowtheretal 



Morris et all (120001) for WN stars and iRepolust et al 



(1995p), 

(l2004l) . lBouret et al. (2005) f or O supergiant s are shown by open 
symbols. The results of lHamann et al.l {2006) (WN stars) are dis- 
played by the black filled symbols. Early WN stars (<WN7) are 
shown by triangles, while late WN stars (>WN8) are circles. 
Ofpe/WN9 stars are star symbols and O supergiants are open 
asterisks. 



jected rotational velocity of the brightest OB supergiants was of 
the ord er 100 km s -1 . Inspection of Fig. 1 of iMevnet & Maederl 
d2003l) reveals that after 4-8 Myrs (age of the GC population of 
early type stars), stars with initial masses in the range 25-60 M 
(appropriate for the OB supergiants) have indeed V sin/ between 
40 and 200 km s _1 (depending on the exact age / mass). This 
is another indication that 1) the evolutionary tracks with rota- 
tion are adequate for the present discussion, and 2) that the GC 
massive stars do not have exceptionally large rotational veloci- 
ties, as could be suspected for such a dense environment where 
interactions should be frequent. 

How does our He determination compare to other analyses 
of massive stars? The answer is given in Fig. [27] Blue circles 
show the He content of Galactic O supergiants. It is important 
to note that the uncertainty usually quoted in such determina- 
tions is ±0.05 - 0.1. The circle with the error bars in Fig. [271 
represents the average OB supergiants in the GC. We see that 
although large, the He content of these stars is still compatible 
with the range of values found for O supergiants. In conclusion, 
the GC OB supergiants are on average similar to other Galactic 
stars in terms of He enrichment. This abundance is in addition 
compatible with present evolutionary tracks. 

iPaumard et al.l (120061) noted that a few OB supergiants had a 
strong He absorption on the blue side of Br y and suggested that 
these stars are significantly He-rich. Unfortunately, given the low 
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Fig. 27. Comparison between derived ratio of He to H abun- 
dance (symbol with error bars) of OB super giants in t he 
Galactic Center to values from other studie s ( Herrer o et 
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S/N ratio spectra of these stars, we could not derive quantitative 
constraints on the He abundance. 



10. Conclusion 

We have carried out a detailed analysis of 18 evolved massive 
stars (Ofpe/WN9 and Wolf-Rayet) located in the central parsec 
of the Galaxy from H and K band spectroscopy obtained with 
SINFONI on the ESO/VLT. The average spectrum of 10 bright 
OB supergiants was also examined. For quantitative analysis, we 
have used state of the art non-LTE atmosphere models including 
winds and line-blanketing computed with the code CMFGEN. 
The main results are: 

• The population of massive stars in the central parsec is able 
to supply the amount of H and Hei ionising photons required 
to reproduce the nebular emission. The contribution of the 
bright, cool post-ma i n sequ enc e massive stars, fi rst discov- 
ered bv lForrest et al.l(ll987l) a nd lAllen et all d 19901) and anal- 
ysed bv lNaiarro et al.l (1 19971) . accounts for only ~4% of the 
total H ionising flux reconciling stellar evolution with obser- 
vations in the Galactic Center. State of the art evolutionary 
and atmosphere models also reconcile the observed stellar 
content with a population synthesis model of a starburst of 
age ~ 6 Myrs. 

• Ofpe/WN9 stars are evolved mass i ve sta rs close to a LBV 
phase. Compared to iNajarro et al.l (1 19971) . we find they are 
less He rich, slightly hotter (although the temperature is 
poorly constrained) and have lower mass loss rates. WN8 
stars are found to have properties similar to other Galactic 
Wolf-Rayet stars of the same spectral type. From morpho- 
logical as well as quantitative arguments, they are likely the 
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descendents of Ofpe/WN9 stars. Two stars are classified as 
WN8/WC9 stars. Their quantitative analysis reveals that, as 
suggested by their spectral type, they show both H and He 
burning products at their surface. A direct evolutionary link 
between the GC Ofpe/WN9, WN8 and WN/C stars of the 
form 

(Ofpe/WN9 ^ LBV) -> WN8 -> WN/C 

is proposed, similar to lCrowther et al.l {l995b). 

• Quantitatively, stellar evolutionary tracks with rotation and 
Z = Z Q overpredict the luminosity of the GC Wolf-Rayet 
stars. Tracks with Z = 2 x Z Q are more appropriate. This 
may indicate that the mass loss rates adopted in the current 
evolutionary tracks during the Wolf-Rayet phase are too low. 
However, accurate metallicity determinations are needed to 
solve this issue. 

• On average, GC OB supergiants are He rich, but not signif- 
icantly richer than other galactic supergiants. Their present 
projected rotational velocities is ~ 100 km s _I . These proper- 
ties are quantitatively compatible with stellar evolution with 
rotation on and close to the main sequence. 

• The wind properties of WN stars in the Galactic center are 
very similar to other Galactic stars of the same spectral type. 
The luminositie s of late WN s t ars ar e lower in the GC than 
in the sample of Hamann et al. ( 2006) but in good agreement 
with the results of lCrowther et al.l dl995bl) . We argue that the 
difference with the results of Hamann et al. is due to uncer- 
tain distances for their sample stars, while in our case the 
distance to the GC is well constrained. 

Our study has shown that the GC massive stars are, on av- 
erage, similar to other Galactic stars. This strongly suggests that 
they follow a common evolution, regardless of their possible dif- 
ferent formation process. A key question which remains to be 
addressed in future studies is the metallicity in the GC. A better 
knowledge of this parameter is crucial to quantitatively test evo- 
lutionary models. It is also crucial in the context of the chemical 
evolution of the Galaxy. Also important is the study of nebular 
emission in the GC by means of recently developed 3D photo- 
ionization models. Our total stellar SED will be a crucial input 
for such models which will likely constrain the geometry and 
density of the ionised gas. 
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